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’An  experimental  and  analytical  study  of  a rotary  vane  steam  expander  was 
conducted  to  determine  the  effect  of  leakage,  friction  and  heat  transfer  on 
the  expander  performance.  A commercially  available  rotary  vane  air  motor 
was  modified  to  operate  on  steam  utilizing  little  or  no  liquid  lubricant. 
The  indicated  power  output,  shaft  power  output  and  frictional  power  loss  of 
the  vane  expander  were  experimentally  determined  as  a function  of  speed, 
inlet  timing  and  supply  steam  conditions.  The  steam  mass  flow  rate  and 
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component  temperatures  were  also  measured.  The  data  show  that  severe 
Internal  leakage  and  frictional  energy  dissipation  were  major  causes 
of  efficiency  reductions. 

An  analytical  model  of  the  expander  thermodynamics,  friction, 
leakage  and  heat  transfer  was  developed  from  fundamental  principles. 
The  model  predicts  the  expander  leakage  flow  rate,  frictional  power 
loss  heat  transfer  rate  and  the  effect  of  these  losses  on  the  power 
output  and  efficiency.  The  analytically  and  experimentally  determined 
frictional  power  losses  were  in  agreement.  The  component  temperature 
profiles  were  predicted  with  maximum  errors  of  10%  - 15%.  The 
predicted  leakage  flow  was  approximately  16%  below  the  experimentally 
determined  value.  Errors  in  the  leakage  flow  predictions  resulted  in 
the  predicted  indicated  power  outputs  being  20%  - 40%  below  the 
experimental  values.  This  was  condiered  reasonably  good  in  light  of 
the  difficulty  in  identifying  the  steam  leakage  paths  and  component 
clearances . 
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ABSTRACT 


An  experimental  and  analytical  8tudy  of  a rotary  vane  steam 
expander  was  conducted  to  determine  the  effect  of  leakage,  friction  and 
heat  transfer  on  the  expander  performance.  A commercially  available 
rotary  vane  air  motor  was  modified  to  operate  on  steam  utilizing  little 
or  no  liquid  lubricant.  The  indicated  power  output,  shaft  power  output 
and  frictional  power  loss  of  the  vane  expander  were  experimentally 
determined  as  a function  of  speed,  inlet  timing  and  supply  steam 
conditions.  The  steam  mass  flow  rate  and  component  temperatures  were 
also  measured.  The  data  show  that  severe  internal  leakage  and  frictional 
energy  dissipation  were  major  causes  of  efficiency  reductions. 

An  analytical  model  of  the  expander  thermodynamics,  friction, 
leakage  and  heat  transfer  was  developed  from  fundamental  principles. 

The  model  predicts  the  expander  leakage  flow  rate,  frictional  power  loss 
heat  transfer  rate  and  the  effect  of  these  losses  on  the  power  output  and 
efficiency.  The  analytically  and  experimentally  determined  frictional 
power  losses  were  in  agreement.  The  component  temperature  profiles  were 
predicted  with  maximum  errors  of  10%  - 15%.  The  predicted  leakage  flow 
was  approximately  16%  below  the  experimentally  determined  value.  Errors 


in  the  leakage  flow  predictions  resulted  in  the  predicted  indicated 
power  outputs  being  20%  - 40%  below  the  experimental  values.  This  was 
considered  reasonably  good  in  light  of  the  difficulty  in  identifying  the 
steam  leakage  paths  and  component  clearances. 
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high  temperature  vapor  cycles  has  been  limited. 


To  determine  the 


applicability  of  vane  expanders  in  Rankine  cycle  applications  and  to 
improve  their  performance,  it  is  necessary  to  understand  the  various 
losses  which  occur.  Consequently,  it  is  the  object  of  this  research  to 
determine  the  effects  of  friction,  leakage,  and  heat  transfer  on  the 
performance  of  a rotary  vane  expander  through  a combined  analytical  and 
experimental  program. 
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1.2  Previous  Related  Studies 

Wolgemuth  and  Olson  [1]  developed  a thermodynamic  model  of  the 
rotary  vane  steam  expander  to  study  its  general  operating  characteristics. 
The  analysis  included  the  transient  charging  and  discharging  processes 
of  the  expander  to  permit  the  breathing  of  these  devices  to  be  studied. 
Although  the  model  contained  provisions  for  the  inclusion  of  heat  transfer, 
and  was  adaptable  to  include  the  effects  of  leakage  and  friction,  a 
detailed  analysis  of  these  phenomena  and  their  effect  on  the  expander 
performance  was  not  performed.  Their  study  shows  that  easy  breathing  in 
vane  expanders  is  readily  achieved  and  that  a vane  expander  operating  in 
the  absence  of  heat  transfer,  leakage  and  friction  can  obtain  a high 
expansion  efficiency. 

The  frictional  forces  existing  between  the  vanes  and  the  stator 
of  a rotary  vane  air  cycle  refrigeration  machine  (ROVAC)  were  analyzed 
by  Edwards  and  McDonald  [2].  The  model  utilized  mean  values  of  the 
geometric  variables  involved  in  the  computation  of  the  frictional  forces. 

In  order  to  determine  the  pressure  forces  on  the  vanes,  the  analysis 
used  a thermodynamic  model  which  did  not  have  the  capability  to  predict 
transient  pressure  effects. 
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Peterson  and  McCahan  [3]  developed  a thermodynamic  model  of  an 
oil  flooded  sliding  vane  air  compressor.  The  model  analyzed  the 
frictional  forces  existing  in  the  compressor  and  used  a dynamic  analysis 
of  the  working  fluid  to  compute  the  pressure  forces  on  the  vanes.  The 
model  included  the  leakage  occurring  between  the  primary  control  volume 
(volume  between  adjacent  vanes)  and  the  control  volumes  leading  it  and 
lagging  it  as  well  as  between  the  primary  control  volume  and  those 
regions  at  inlet  pressure.  In  computing  the  leakage  flow,  an  empirically 
determined  discharge  coefficient  was  employed  to  account  for  the  blockage 
of  the  flow  area  by  oil.  The  heat  transfer  was  modeled  by  assuming  the 
air  transferred  heat  only  with  the  oil  in  the  chamber.  Good  agreement 
was  obtained  between  the  experimental  and  computed  power  input  require- 
ments and  air  flow  rates. 

Eckard  [4]  graphically  presented  experimental  data  taken  from  a 
high  temperature  vane  expander  being  developed  by  General  Electric. 

Eckard  pointed  out  that  leakage  and  friction  are  two  major  problem  areas 
encountered  in  the  development  of  a multi-vane  vapor  expander,  with 
practical  solutions  to  leakage  problems  the  most  difficult  to  obtain. 

In  order  to  reduce  leakage  without  large  increases  in  friction,  the 
General  Electric  expander  depends  on  the  maintenance  of  small  clearances. 
Eckard  noted,  however,  that  considerable  care  must  be  taken  to  achieve 
and  maintain  small  clearances  under  the  changing  temperature  conditions 
encountered.  Additional  observations  made  by  Eckard,  based  on  the  data 


collected  are: 
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1)  Moderate  shift  in  expander  speed  does  not  severely 
reduce  expander  efficiency. 

2)  Breathing  losses  in  a vane  expander  can  essentially 
be  eliminated  since  dynamic  valves  are  not  required. 

3)  Low  speed  operation  (1200-1800  RPM)  is  currently 
essential  for  low  frictional  power  losses. 

Eckard  concluded  that  additional  material  and  design  innovations  are 
necessary  to  further  improve  the  expander  performance  and  to  increase 
its  applicability. 

1.3  Purpose 

The  preceedlng  discussion  indicates  that  vane  expanders  have  a 
potential  for  high  efficiency  in  vapor  cycle  applications  provided 
degradation  of  expander  performance  due  to  friction,  heat  transfer  and 
leakage  can  be  controlled. 

Some  experimental  work  has  been  done  to  develop  an  efficient 
rotary  vane  expander.  Additionally,  friction  and  to  a lesser  extent 
leakage,  have  been  modeled  in  vane  air  compressors  and  in  the  Rovac 
machine  utilizing  techniques  applicable  to  a vane  expander.  However,  a 
physical  description  of  the  combined  phenomena  of  heat  transfer,  leakage 
and  friction  that  occur  within  a multi-vane  vapor  expander  has  not  been 
formulated . 

Since  the  amounts  of  leakage,  heat  transfer  and  a friction  are 


not  independent  of  one  another,  and  since  it  is  difficult,  at  times, 
to  determine  which  of  these  phenomenon  is  associated  with  an  observed 
degradation  of  power  or  efficiency,  an  understanding  and  description  of 
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the  physical  processes  occurring  in  the  expander  is  required  to  identify 
the  mechanisms  of  the  various  losses.  Furthermore,  the  development  of  a 
working  model  of  the  vane  expander  would  permit  a degree  of  optimization 
of  the  expander  design  to  be  achieved  prior  to  the  construction  of  hard- 
ware. 

In  view  of  the  high  efficiency  potential  of  vane  expanders,  and 
in  view  of  the  necessity  to  obtain  a better  understanding  of  the 
phenomena  occurring  within  them,  the  following  are  the  specific 
objectives  of  this  thesis. 

1.  Analytical 

The  objective  of  the  analytical  study  is  to  develop  a working 
model  of  multi-vane  expanders  to  permit  the  computation  of  the  following, 
for  various  fluid  conditions  and  expander  geometry: 

a.  Frictional  power  loss. 

b.  Leakage  flow  rates  and  the  associated  effects  on 
power  output  and  efficiency. 

c.  Heat  transfer  rates  and  their  effects  on  power 
output  and  efficiency. 

d.  Expander  power  and  mass  flow  rate  characteristics 
in  the  ideal  case  and  under  the  influence  of 
leakage,  heat  transfer  and  friction. 

2.  Experimental 

The  objectives  of  the  experimental  work  are  as  follows: 

a.  Measure  the  mass  flow  rate,  and  power  characteristics 
of  a multi-vane  expander  as  a function  of  speed, 
and  arc  of  admission  to  compare  with  the  predicted 
values. 
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b.  Determine  the  frictional  power  loss  as  a function 


of  speed,  and  arc  of  admission  to  obtain  a better 


understanding  of  the  friction  in  vane  expanders 


and  to  test  the  validity  of  the  friction  analysis. 


c.  Determine  the  overall  leakage  flow  rate  in  order 


to  verify  the  leakage  analysis. 


d.  Measure  the  rotor  and  stator  temperatures  for 


developing  and  verifying  the  heat  transfer 


analysis  and  to  determine  the  approximate  amount 


of  rotor  and  stator  thermal  expansion. 
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CHAPTER  II 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 
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2 . 1 Description  of  Expander 

The  vane  expander  used  to  obtain  the  experimental  data  is 
schematically  illustrated  in  Figure  1.  The  configuration  of  the  expander 
was  the  result  of  extensive  modification  to  a commercially  available  Gast 
Corporation  model  8AM  air  motor.  The  expander  primarilly  consisted  of  a 
cylindrical  rotor,  containing  eight  sliding  vanes,  eccentrically  mounted 
in  a cylindrical  housing.  Figure  2 is  a photograph  showing  an  exploded 
view  of  the  expander.  To  permit  the  expander  to  operate  with  little  or 
no  liquid  lubricant  which  might  contaminate  the  working  fluid,  the  vanes 
and  end  plates,  shown  in  Figure  2,  were  constructed  of  Pure  Carbon  P5N 
and  P9  carbon,  respectively.  To  provide  a hard  rubbing  surface  for  the 
vanes,  the  expander  housing  contained  a heat  treated  416  stainless-steel 
liner,  located  in  the  housing  by  two  keys.  The  inlet  and  exhaust  port 
timing  was  controlled  by  the  location  of  the  ports  in  the  liner.  The 
exhaust  port  was  positioned  so  that  exhaust  began  when  the  volume  between 
adjacent  vanes  (the  control  volume)  was  a maximum.  A 135°  arc  of  exhaust 
was  employed  so  that  the  control  volume  exhausted  from  its  maximum  volume 
to  its  minimum  volume  position.  Ideally,  this  eliminated  recompression 
losses. 
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Admission  of  the  steam  was  Initiated  when  the  volume  between 
adjacent  vanes  was  a minimum  (immediately  after  the  exhaust  port  was 
closed).  Two  arcs  of  admission,  22.5°  and  45°,  were  utilized  in  the 
test  program.  Initially,  a 22.5°  arc  of  admission  was  employed  since 
it  permitted  a higher  expansion  ratio  to  be  obtained  than  the  45°  arc 
of  admission.  The  expansion  ratio  is  defined  as  the  ratio  of  the  volume 
at  exhaust  port  opening  to  the  volume  at  intake  port  closing.  For  a 
22.5°  arc  of  admission.  Figure  3 shows  that  in  the  absence  of  friction, 
leakage,  and  heat  transfer,  the  expansion  of  the  working  fluid  from  the 
supply  pressure,  of  150  psia,  to  the  exhaust  manifold  pressure 
(atmospheric  pressure)  is  practically  complete  prior  to  the  exhaust  port 
opening.  Ideally,  this  results  in  a high  adiabatic  expansion  efficiency. 
The  adiabatic  expansion  efficiency  is  defined  as  the  ratio  of  the  work 
per  pound  of  fluid  to  the  isentropic  enthalpy  change  from  the  inlet 
manifold  conditions  to  the  exhaust  manifold  pressure.  In  the  ideal  case, 
the  45°  arc  of  admission  results  in  an  increase  in  the  expander  power 
output,  as  indicated  by  the  increase  in  area  under  the  P*  vs.  V*  diagram 
shown  in  Figure  3.  However,  the  increased  arc  of  admission  also  results 
in  a decrease  in  the  expansion  ratio  and  hence  a decrease  in  adiabatic 
expansion  efficiency. 

To  control  leakage  between  the  ends  of  the  rotor  and  the  end 
plates,  shims  were  installed  between  the  carbon  end  plates  and  the 
housing.  The  shims  permitted  the  cold  clearance  between  the  rotor  and 
end  plate  to  be  adjusted  to  obtain  a balance  between  leakage  and  end 
plate-rotor  friction. 
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SUPPLY  PRESSURE  I50PSIA 
SPEED  1007  RPM 

• 22.5*  ARC  OF  ADMISSION 


0.0 

0.0 


l I I I I 

0.05  0 

V*  (DIMENSIONLESS  VOLUME) 


Figure  3.  Dimensionless  Pressure  Versus  Dimensionless  Volume, 

Ideal  Expander  (No  Leakage,  Heat  Transfer,  or  Friction) 
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Contact  between  the  vanes  and  the  stator  was  maintained  by  a 
combination  of  springs  and  steam  pressure  acting  at  the  base  of  the 
vanes.  Steam  was  ported  to  the  base  of  the  vanes  by  two  1/8  inch  and 
one  3/16  inch  diameter  holes  (pressure  balancing  ports  or  pressurization 
ports)  which  extended  from  the  control  volume  leading  a vane  to  the  vane 
base,  as  shown  schematically  in  Figure  1.  The  spring  force  on  the  vane 
base  was  maintained  by  a combination  of  push  pins  and  leaf  springs.  The 
springs  had  a spring  constant  of  approximately  729  lb^/in.  The  variation 
in  the  force  exerted  by  the  spring  on  the  vane  is  shown  in  Figure  4.  A 
high  spring  force  was  maintained  on  the  vane  when  the  vane  was  near  the 
inlet  port.  This  was  necessary  to  help  prevent  the  high  pressure  steam 
from  pushing  the  vane  into  the  vane  slot.  Lower  spring  forces  were 
exerted  on  the  vane  during  portions  of  the  closed  expansion  and  the 
exhaust  processes  (where  the  steam  pressures  were  lower)  to  reduce 
frictional  energy  dissipation.  To  develop  a large  pressure  force  on  the 
vane  base  as  it  moved  from  the  exhaust  side  to  the  inlet  side  of  the 
expander,  prepressurization  ports  were  located  in  the  end  plates  as 
shown  in  Figure  1.  The  ports  permitted  the  flow  of  high  pressure  steam 
to  the  vane  base  while  the  vane  tip  was  exposed  to  the  low  exhaust 
pressure.  This  method  of  obtaining  a vane-stator  seal  resulted  from 
testing  the  foll<^ing  sealing  configurations: 

1)  Spring  force  only  acting  on  the  vane  base  (Spring  constant 
of  17.1  lb^/in.  and  729  lb^/in.).  In  both  cases,  the  vane 
did  not  seal  against  the  outer  stator. 

2)  Steam  pressure  only,  acting  on  the  vane  base  utilizing 
pressure  balancing  ports.  Pressurization  of  the  vane  base 
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at  the  time  of  inlet  opening  was  not  rapid  enough  to  prevent 
the  vane  from  being  pushed  into  the  slot  by  high  pressure 
supply  steam.  This  resulted  In  severe  leakage. 

3)  Steam  pressure  only  acting  on  the  vane  base  utilizing  pressure 
balancing  and  prepressurization  ports.  The  prepressuriza- 
tion ports  were  not  effective  in  rapidly  pressurizing  the 
base  of  the  vane.  This  also  resulted  in  severe  leakage  since 
the  vane  was  pushed  into  the  slot  by  the  high  pressure  steam. 

4)  Combined  spring  force,  with  a spring  constant  of  729  lbf/in., 
and  steam  pressure  utilizing  pressure  balancing  and  pre- 
pressurization ports  acting  on  the  vane  base. 


2.2  Experimental  Apparatus 


A schematic  diagram  of  the  apparatus  used  to  test  the  expander  is 
shown  in  Figure  5.  The  expander  was  coupled  to  a General-Electric  D-C 
cradled  dynamometer  model  number  5B284B1010,  to  obtain  torque  and  power 
measurements.  The  dynamometer  force-measuring  device  was  a Toledo 
balance  with  a 0-30  lb^  range  in  0.01  lbm  subdivisions.  The  motor  speed 
was  measured  with  a General-Electric  model  An5531-1  tachometer  generator 
coupled  to  a Standard  Electric  Time  Company  type-SG6  RPM  counter  and  a 
clock. 


A throttling  valve  in  the  steam  supply  line  permitted  the  supply 
steam  to  be  throttled  from  250  psia  to  150  psia  or  less  before  entering 
the  expander.  To  ensure  superheating  of  the  supply  steam,  three 
Cole-Palmer  flexible  heating  tapes,  capable  of  delivering  2.88  kW  at 
230  volts  were  wrapped  around  the  steam  supply  line.  Control  of  the 
energy  input  into  the  heating  tapes  was  obtained  by  use  of  a Variac. 


Figure  5.  Schematic  of  Experimental  Apparatus 
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I Omega  Engineering  Company  chromel-alumel  thermocouples  were 

I mounted  in  the  steam  lines  to  measure  the  inlet  and  exhaust  temperatures. 

A mixing  tank  was  installed  in  the  inlet  line  to  eliminate  radial 
temperature  gradients  in  the  flow  prior  to  making  the  inlet  temperature 
measurements.  This  permitted  a better  determination  of  the  inlet  mixing 

(cup  temperature.  The  inlet  and  exhaust  temperatures  were  continuously 

recorded  on  a Leeds  and  Northrup  Type  H potentiometric  recorder. 

The  expander  inlet  and  exhaust  pressures  were  measured  by  ( 

t 

Ashcrof t-Bourdon  tube  pressure  gages,  0-200  psig  and  0-100  psig,  i 

respectively.  1 

The  steam  mass  flow  rate  was  measured  by  condensing  the  steam  in  j 

? 

a tank  of  water  and  measuring  the  Increase  in  tank  weight  over  a known  ! 

I ; ; 

j time  interval.  The  combined  weight  of  the  tank  and  the  water  was  j i 

measured  with  a 0-1500  lb  Toledo  scale  having  0.5  lb  subdivisions. 

ra  m 

(The  tare  weight  of  the  tank  was  measured  with  the  tank  empty  so  that 
| the  weight  of  the  water  could  be  obtained).  A three  way  valve  was 

■ installed  between  the  expander  outlet  and  the  tank  inlet  so  that  the 

i 

tank  could  be  bypassed  until  the  mass  flow  rate  measurements  were  made. 

Since  the  inlet  steam  was  superheated,  its  state  was  determined 
from  the  pressure  and  temperature  measurements.  However,  since  the 
exhaust  steam  could  be  a two  phase  mixture,  provisions  were  made  for 
determining  the  enthalpy  of  the  exhausting  steam  by  making  an  energy 
balance  on  the  condensing  tank  used  to  measure  the  mass  flow  rate.  To 
do  this,  the  tank  was  heavily  insulated.  The  tank  water  temperature  was 
measured  using  a thermistor  (Yellow  Springs  Instrument  Company  No.  701) 
with  a digital  readout  (United  Systems  Corporation  Model  581C).  The 

I 
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system  accuracy  was  jK).55°C.  Additionally,  chromel-alumel  thermocouples 
were  submerged  in  the  water  at  different  depths  to  provide  a check  on 
temperature  gradients  in  the  water.  The  tank  wall  temperature  was 
measured  with  five  (5)  chromel-alumel  thermocouples  mounted  at  various 
vertical  positions  on  the  tank  wall.  The  output  of  the  thermocouples 
was  measured  with  a Leeds  and  Northrup  model  8686  potentiometer. 
Measurement  of  the  steam  flow  into  the  tank,  water  temperature  rise,  and 
tank  wall  temperature  rise  permitted  an  energy  balance  to  be  performed 
on  the  tank.  The  enthalpy  of  the  two  phase  exhausting  steam  could  then 
be  computed  from  the  following  expression: 
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It  was  discovered  during  expander  tests  that  leakage  in  the 
expander  was  sufficiently  high  so  that  the  exhausting  steam  was  super- 
heated. Therefore,  its  state  was  determined  by  the  exhaust  temperature 
and  pressure  measurements. 

Oil  was  pumped  to  the  expander  by  a Mandel  class  XN  force  feed 

lubricator.  The  lubricator  was  driven  by  a variable  speed  motor.  The 

oil  flow  rate  was  controlled  by  the  speed  of  the  motor  and  the  stroke 

length  of  the  pump  metering  plunger.  The  pump  was  capable  of  delivering 

oil  flows  of  from  0 1b  /min.  to  0.13  lb  /min.  The  oil  employed  was  mixed 

m m 


using  the  raw  materials  and  quantities  given  in  Reference  3. 
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Measurement  of  the  control  volume  pressure  was  accomplished  with 
three  piezoelectric  pressure  transducers  mounted  in  one  end  plate.  The 
positions  of  the  transducers  permitted  measurement  of  the  control  volume 
pressure  during  filling  of  the  volume,  during  expansion  and  during  the 
initial  phase  of  control  volume  blowdown.  Two  Kistler  model  603  pressure 
transducers  coupled  to  Kistler  model  504  charge  amplifiers  were  used. 

The  third  transducer  was  a Metrix  model  5016  connected  to  a Metrix  model 
5080-3  charge  amplifier.  Provisions  were  made  for  simultaneously 
monitoring  the  charge  amplifier  outputs  with  a Tektronix  model  5103N 
oscilloscope  and  a CEC  (Consolidated  Electrodynamics  Corporation)  model 
5-124A  oscillograph. 

Since  the  pressure  transducers  were  not  designed  to  withstand  the 
high  temperature  steam  environment  and  since  their  output  was  temperature 
sensitive,  the  transducers  were  connected  to  water  cooling  adapters.  The 
adapters  permitted  the  output  of  the  transducers  to  remain  stable. 

Seven  (7)  chrome 1- a lumel  thermocouples  (0.005  inches  in  diameter) 
were  mounted  around  the  periphery  of  the  stainless  steel  liner.  The 
thermocouples  were  positioned  0.050  inches  from  the  inside  surface  of 
the  liner.  The  thermocouple  leads  were  taken  out  of  the  expander 
housing  through  conax  glands. 

Chromel-alumel  thermocouples  were  also  installed  at  eleven 
circumferential  positions  in  the  carbon  end  plate.  The  thermocouples 
were  mounted  0.065  inches  from  the  inner  face  of  the  end  plate.  The 
liner  and  end  plate  thermocouples  were  connected  through  a switch  and  a 
CEC  Type  1-165  amplifier  to  the  oscillograph. 
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The  surface  temperature  of  the  rotor  was  measured  by  an  iron- 
Conscantan,  Bendersky  [6]  type  thermocouple  having  a response  time  of 
1 microsecond.  The  thermocouple  is  shown  mounted  in  Figure  6.  The 
signal  from  the  thermocouple  was  transferred  from  the  rotor  through  a 
mercury  slip  ring  (Meridian  Laboratory  Incorporated,  Mercotac  I).  The 
signal  was  amplified  by  a CEC  Type  1-165  high-gain  amplifier  and 
recorded  on  the  oscillograph.  The  Mercury  slip  ring  was  capable  of 
operating  at  a maximum  temperature  of  150°F.  To  ensure  that  the  slip 
ring  temperature  did  not  exceed  this  value,  a cooling  adapter  cooled  an 
extension  of  the  rotor  shaft,  ahead  of  the  slip  ring.  This  permitted 
the  slip  ring  to  operate  at  room  temperature. 

2 . 3 Experimental  Procedure 

Prior  to  operating  the  expander,  the  steam  line  was  heated  by  the 
heating  tapes.  The  steam  flow  was  diverted,  ahead  of  the  expander,  to 
a drain  line  to  remove  any  liquid  condensate  in  the  steam  line  (see 
Figure  5).  Once  the  condensate  was  removed,  the  solenoid  valve  upstream 
of  the  expander  was  opened,  causing  the  expander  to  rotate,  and  the 
drain  valve  was  closed.  The  force  feed  lubricator  was  started,  and 
cooling  water  flow  to  the  pressure  transducer  and  slip  ring  cooling 
adapters  was  initiated.  The  supply  pressure  was  adjusted  to  the  desired 
value,  utilizing  the  throttling  valve.  The  dynamometer  load  was  adjusted 
until  the  expander  was  rotating  at  the  desired  speed.  The  steam 
exhausting  from  the  expander  bypassed  the  condensing  tank  via  a three 
way  valve. 

The  expander  was  operated  in  this  manner  until  the  inlet  steam 
temperature  and  the  expander  component  temperatures  reached  steady  state 
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Figure  6,  Bendersky  Thermocouple  Mounting 
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values.  The  expander  was  then  shut  down  by  closing  the  solenoid  and  the 
throttling  valves.  The  force  feed  lubricator  was  also  stopped.  The 
initial  mass  of  water  in  the  condensing  tank  was  recorded  and  the 
pressure  transducer  outputs  grounded.  The  solenoid  valve  was  opened  and 
the  throttling  valve  adjusted  to  provide  the  correct  supply  pressure. 

The  force  feed  lubricator  was  started  and  the  revolution  counter  and 
clock  were  activated.  The  expander  exhaust  steam  was  collected  in  the 
condensing  tank.  The  thermocouple  and  pressure  transducer  outputs  were 
recorded  along  with  the  supply  and  exhaust  pressures.  The  dynamometer 
load  was  also  recorded.  After  approximately  three  to  five  minutes  of 
operation,  the  revolution  counter  and  clock  were  stopped,  the  condensing 
tank  was  bypassed,  the  force  feed  lubricator  was  stopped  and  the  solenoid 
and  throttling  valves  were  closed.  The  final  mass  of  water  in  the 
condensing  tank  was  recorded  along  with  the  number  of  revolutions  and 
the  operating  time.  Water  flow  to  the  cooling  adapters  was  continued 
until  the  next  test  or  until  the  expander  components  had  cooled  to  room 
temperature. 

2.4  Error  Analysis 

The  errors  associated  with  the  various  measurements  have  been 
estimated  by  considering  the  parameters  which  affect  the  absolute 
accuracy  of  the  measurements. 


2.4.1  Power  Measurement . Errors  in  the  power  measurements 
resulted  from  an  uncertainty  in  the  length  of  the  torque  arm,  the  motor 
speed  and  the  torque  arm  loading.  The  estimated  errors  associated  with 


the  power  measurements  are  tabulated  in  Table  1. 
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An  overall  estimate  of  the  error  in  the  power  measurement  was 
made  by  taking  the  square  root  of  the  sum  of  the  squares  of  the  errors 
given  in  Table  1.  The  resulting  estimated  error  is  approximately  +0-01  hp. 

2.4.2  Mass  Flow  Rate  Measurement.  Calibration  of  the  Toledo 

scale  showed  it  to  be  accurate  to  within  0.5  lb  . An  error  of  0.5  lb 

m m 

in  measuring  the  mass  of  water  in  the  tank  resulted  in  an  error  in  the 
mass  flow  rate  measurement  of  approximately  +0.1  lb^/min. 

2.4.3  Control  Volume  Pressure  Measurement.  An  estimate  of  the 
error  in  the  control  volume  pressure  measurement  is  quite  difficult. 

The  output  of  the  piezoelectric  transducers  is  affected  by  transducer 
temperature  and  decay  of  the  output  signal.  The  control  volume  pressure, 
measured  by  the  pressure  transducers  during  the  portion  of  the  cycle 
when  the  inlet  was  open,  was  compared  to  the  supply  pressure  gage 
measurements.  General  agreement  was  within  +4  psia.  Similarly,  agree- 
ment between  overlapping  pressure  transducers  which  simultaneously 
measured  the  same  control  volume  pressures  during  certain  portions  of 
the  cycle  was  within  approximately  +4  psia.  Hence,  a rough  estimate  of 
the  error  associated  with  the  pressure  measurement  is  +4  psia. 

2.4.4  Steam  Enthalpy  Measurement.  The  errors  associated  with 
the  measurement  of  the  expander  inlet  and  exhaust  steam  enthalpy  result 
from  errors  in  the  steam  temperature  and  pressure  measurements.  Errors 
also  result  from  stratification  of  the  flow.  A mixing  tank  was  installed 
in  the  inlet  steam  line  to  minimize  flow  stratification.  The  superheated 
exhaust  steam  was  mixed  during  the  control  volume  blowdown  process.  The 
exhaust  temperature  and  pressure  measurements  were  made  adjacent  to  the 
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expander  exhaust  port  so  that  stratification  of  the  flow  in  the  steam 
i line  would  have  a very  short  time  to  occur. 

Errors  in  the  pressure  measurement  of  +2  psia  and  in  the  tempera- 
ture measurement  of  +4°F  result  in  errors  in  the  steam  enthalpy  of  +3.19 

BTU/lb  . This  is  approximately  a +0.3%  error.  Therefore,  in  view  of 
m 

t he  precautions  taken  to  eliminate  flow  stratification,  an  error  of 
' +1%  in  the  enthalpy  measurement  can  be  expected  with  reasonable  confidence. 

2.4.5  Temperature  Measurement.  The  manufacturer  specified 
deviation  of  the  chromel-alumel  thermocouples,  output-vs-temperature 
characteristics  from  the  standard  table  values  is  +4°F.  This  was 
verified  through  experimental  tests. 

Calibration  tests  of  the  iron-Constantan , Bendersky  thermocouple 

i 

showed  the  deviation  of  its  output-vs-temperature  characteristics,  from 
the  standard  table  values  to  be  +4°F. 

Precautions  were  taken  to  reduce  thermocouple  installation  errors. 

) Particular  care  was  taken  to  minimize  the  change  in  the  rotor  heat 

, transfer  characteristics  caused  by  the  installation  of  the  Bendersky 

surface  thermocouple.  An  iron-Constantan  thermocouple  was  employed. 

The  thermocouple  had  an  iron  body  so  that  it  had  thermal  properties 
similar  to  those  of  the  iron  rotor.  Thin  thermocouple  wires  were  used 
to  minimize  lead  losses.  The  thermocouple  was  mounted  so  that  the 
pattern  of  the  steam  flow  over  the  surface  was  disturbed  as  little  as 

possible  by  the  thermocouple. 

it 


It  is  believed  that  the  precautions  taken  in  mounting  the 

! 

[ Bendersky  thermocouple  reduced  the  thermocouple  installation  errors. 


CHAPTER  III 
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ANALYTICAL  CONSIDERATIONS 
3.1  Thermodynamic  Model 

The  foundation  of  the  vane  expander  analysis  is  the  thermodynamic 
model  described  by  Wolgemuth  and  Olson  [1].  The  thermodynamic  analysis 
utilized  the  principles  of  conservation  of  energy  and  mass,  as  applied 
to  the  volume  between  adjacent  vanes  (the  control  volume),  in  combination 
with  the  geometric  expressions  for  the  volume,  and  the  equation  of  state 
in  the  form  of  steam  tables.  These  relationships  were  used  to  obtain 
differential  equations  for  pressure,  mass,  volume  and  in  the  superheated 
region,  temperature  of  the  working  fluid  as  a function  of  rotor  position. 
The  differential  equations,  in  conjunction  with  the  expression  for  the 
indicated  power  output  of  the  control  volume,  10  * PV  , were  solved 
numerically  using  Hamming's  Predictor  Corrector  method  [7J. 

In  formulating  the  thermodynamic  model,  the  first  law  of  thermo- 
dynamics was  employed  in  the  form 

Q - W + (mh) . - (mh)  = U . (3-1> 

in  out  cv 

Therefore,  the  model  is  capable  of  predicting  the  power  output  of  the 
expander  in  the  ideal  case,  where  mass  transfer  occurs  only  through  the 
inlet  and  the  exhaust  ports,  and  the  control  volume  is  adiabatic. 

However,  the  thermodynamic  model  can  also  determine  the  effects  of  heat 
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transfer  on  the  expander  performance  by  computing  the  heat  transfer  rates 

in  separate  heat  transfer  models  and  inputting  the  rates  into  the 

thermodynamic  model  as  Q . Similarly,  the  effect  of  leakage  can  be 

determined  by  computing  the  leakage  flow  in  a separate  leakage  model  and 

inputting  the  leakage  flow  rates  into  the  thermodynamic  model  as  m 

or  m ^ . Hence,  in  the  case  of  leakage,  m,  or  m includes  flow 
out  0 in  out 

through  the  inlet  and  exhaust  ports  as  well  as  leakage. 

The  work  term  in  Equation  (3.1)  is  the  indicated  work  of  the 
expander  or  the  PdV  work  done  by  the  working  fluid.  The  effect  of  fric- 
tion on  the  shaft  power  can  be  determined  by  computing  the  frictional  < 

work  loss  by  a separate  friction  model  and  subtracting  it  from  the 
indicated  work. 

To  analytically  determine  the  characteristics  of  the  expander  in 
the  presence  of  leakage,  friction,  and  heat  transfer,  separate  models 
for  these  phenomena  were  developed  which  interact  with  the  thermodynamic 
model  as  described  earlier  and  as  schematically  illustrated  in  Figure  7. 

3.2  Friction  Model 

The  analysis  of  the  expander  power  loss  due  to  friction  assumes 
that  there  is  contact  and  hence,  friction  between  the  vanes  and  the 
following  components: 

a.  Stator 

b.  Rotor 

c.  End  Plates. 

There  may  also  be  friction  between  the  rotor  and  the  end  plates 
when  rotor-end  plate  contact  occurs.  It  was  assumed  that  there  was  no 
rotor-end  plate  contact  in  the  experimental  expander.  This  assumption 
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was  based  on  the  observed  freedom  of  rotor  motion  while  the  expander 
was  both  hot  and  cold.  Thus,  it  was  not  included  in  the  analysis. 
However,  since  a small  amount  of  liquid  lubricant  was  employed,  an 
analysis  similar  to  that  described  in  Reference  8 was  used  to  determine 
the  power  loss  resulting  from  the  viscous  forces  acting  on  the  ends  of 
the  rotor.  By  neglecting  centrifugal  effects  and  assuming  the  flow  to 
be  steady,  laminar  and  incompressible,  the  following  expression  for  the 
power  loss  due  to  these  viscous  forces  was  obtained: 

2 4 3 

. UOJ  it  R.  2u)E  R. 

U = — — + — — H£.  (3  2) 

VISCOUS  C 3 dX  • ^ ' 


Reference  9 states  that  the  pressure  drop  term  in  Equation  (3.2)  is 
usually  negligible  since  the  leakage  is  often  normal  to  the  velocity  in 
the  regions  where  viscous  dissipation  occurs.  Therefore,  that  term  was 
neglected. 

In  order  to  determine  the  frictional  forces  on  a vane,  free  body 
diagrams  like  that  shown  in  Figure  8 were  employed. 

The  forces  that  are  the  source  of  the  frictional  energy  dissipa- 
tion are  F , F , F_  , and  F_  . Prior  to  determining  the 
Kd  Rl  3 £ 


magnitude  of  these  forces,  it  was  necessary  to  determine  the  pressure 
forces  and  the  normal  and  tangential  accelerations. 

The  normal  and  tangential  accelerations  of  the  vanes  may  be 
computed  from  the  expander  geometry  under  the  assumption  that  the  vane 
tip  is  in  contact  with  the  stator  at  all  times.  The  expressions  used 
in  computing  the  accelerations  are  given  in  Appendix  B. 


Figure  8.  Vane  Free  Body  Diagram 


The  net  pressure  force  on  the  side  of  the  vane,  , is  computed 


as  follows 


(PLAG  " PLEAD)  A 


(3.3) 


where  P^G  and  Px.EAD  are  t*ie  Pressures  control  volumes  lagging 

and  leading  the  vane  and  A is  the  exposed  area  of  the  vane.  The 
pressures  used  are  those  predicted  by  the  thermodynamic  model  in  the 
presence  of  leakage  but  in  the  absence  of  friction.  Due  to  the  small 
vane-rotor  clearance  in  a cold  expander,  it  was  assumed  that  acts 

only  on  the  vane  in  the  space  between  the  rotor  and  the  stator  and  no 
other  pressure  forces  act  on  the  side  of  the  vane.  In  computing  the 
force  at  the  base  of  the  vane,  it  was  assumed  that  the  pressure  under 
the  vane  was  the  same  as  that  in  the  leading  control  volume  due  to  the 
existence  of  the  pressure  balancing  ports.  The  spring  force  exerted  on 
the  vane  is  equal  to  the  product  of  the  spring  constant  and  the  spring 
deflection. 


Due  to  its  larger  coefficient  of  thermal  expansion,  the  rotor  was 

assumed  to  be  longer  than  the  vanes.  Therefore,  the  force  exerted  on 

the  vanes  by  the  end  plates  was  computed  by  assuming  that  the  vane  was 

pushed  against  one  end  plate  by  the  pressure  in  the  control  volume.  The 

end  force  on  the  vane  was  then  equal  to  the  product  of  the  control  volume 

pressure  and  the  area  of  the  vane  end. 

Since  the  mass  of  the  vanes  is  known,  and  since  the  friction 

coefficient  may  be  estimated,  the  only  unknowns  are  F _ , F _ and  F . 

RT  KB  J 


These  are  found  by  applying  the  three  equations  of  dynamic  equilibrium 


Efnormal 

Vn  • 

(3.4) 

Zftangential 

* n^t  * and 

(3.5) 

Ema  - 0 * 

(3.6) 

to  the  free  body  diagrams  of  the  vane. 

During  one-half  of  the  rotor  revolution,  the  vane  is  pushed  inward 
(toward  the  center  of  the  rotor)  by  the  stator.  During  the  other  half, 
the  vane  is  pushed  outward  by  the  pressure  and  spring  forces  at  its  base. 
Hence,  at  the  end  of  one-half  revolution,  there  is  a reversal  of  the 
direction  of  the  friction  forces  between  the  vane  and  rotor.  Addition- 
ally, the  net  pressure  force  on  the  side  of  the  vane  may  not  always  act 
in  the  direction  of  rotation  as  shown  in  Figure  8.  Therefore,  four  free 
body  diagrams  have  been  used,  resulting  in  four  sets  of  equations.  As 
an  example,  the  equations  for  Figure  8 in  matrix  form  are 


-1  -u2  -p2 

F3 

-PAAEVN-SPF+yiFE-MvAN+WN+PsAP 

-u3  l -1 

X 

frb 

= 

-WE-VVt 

TV  TV 

-lvu3  - — U2  —u2-xv 

frt 

u -p  /LVv  _ __  HRv 

'wt  (T)+MiFE(T)-pf (LV_  +¥t  (T} 

(3.7) 


The  remaining  equations  are  similar  and  will  not  be  presented. 

The  solution  of  the  equations  for  F^  , FRfi  and  FRT  is 
accomplished  by  using  the  method  of  Gaussian  elimination. 
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i The  values  of  F , F , F.,^  , and  F_  and  hence,  the 

I corresponding  friction  forces  were  determined  as  a function  of  0 as  the 

rotor  rotated  in  small  increments  A0  . The  incremental  frictional  work 
i loss  was  computed  as  the  product  of  the  frictional  force  times  the  radial 

or  tangential  distance  the  vane  traveled  during  the  increment  of  rotation. 
The  sum  of  the  incremental  work  losses  over  360°  is  the  total  frictional 
work  loss  per  revolution  due  to  vane  friction.  The  instantaneous  fric- 
tional power  loss  was  computed  at  each  vane  position  as  the  product  of 
the  frictional  force  at  the  position  and  the  velocity. 

The  power  loss  due  to  friction  between  the  vanes  and  rotor,  stator, 
and  end  plate  is  equal  to  the  product  of  the  number  of  vanes,  the 
frictional  work  loss  per  revolution  per  vane  and  the  speed.  The  net 

t 

frictional  power  loss  is  the  sum  of  the  power  loss  due  to  the  vane 
friction  and  the  power  loss  due  to  viscous  drag  on  the  rotor  ends. 

It  was  assumed  that  all  the  energy  dissipated  by  friction  is 
| converted  to  heat  which  is  transferred  into  either  the  rotor,  lubricating 

I 

oil,  stator,  vane  or  end  plate.  A further  discussion  of  the  frictional 

' I 

heat  generation  will  be  given  in  the  description  of  the  heat  transfer 
model. 

3. 3 Leakage  Model 

The  leakage  flow  paths,  assumed  for  this  analysis,  are  shown  for 
a single  control  volume  in  Figure  9. 

It  was  assumed  that  the  leakage  flow  would  be  quasi-steady  and, 
due  to  the  small  clearances  maintained  in  a vane  expander,  laminar. 
Considerable  simplification  was  made  by  treating  the  steam  as  an 
incompressible  fluid.  The  relationship  for  parallel  flow  between  flat 
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plates  was  employed  to  obtain  the  following  expression  for  the  leakage 
flow  along  each  path: 


\ 
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2(P1  ~ P2)b3w 
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where 


3AXp 

,,  3 
2b  w 


(3.8) 


The  lengths  of  the  leakage  paths,  Ax  in  Equation  (3.8)  were 
taken  as  the  straight  line  distance  from  the  source  to  the  sink  along 
the  leakage  path.  The  widths  of  the  leakage  paths  (w)  were  determined 
by  apportioning  the  total  width  of  a leakage  path  at  the  source  over  the 
total  number  of  sinks  the  source  leaked  to  along  the  leakage  path.  It 
should  be  noted  that  Equation  (3.8)  neglects  relative  motion  between  the 
moving  and  stationary  parts. 

Equation  (3.8)  was  applied  to  all  the  leakage  paths  originating 
from  each  control  volume.  As  an  example,  for  control  volume  5 (see 
Figure  9),  the  expression  for  the  net  leakage  flow  out  is 


\5 


(P5  - V<E  r7T>  + <P5  - V<£  rTT>  + <p5  - V iTT 
5-4  5-14  5-3 


+ (P5  ' P2>  + (P5  " Pl)  VT  + <P5  - V + lP5  - V 5^7 


+ <P5  - V<£  + (P5  - P12>  V77  + "5  - PH>  VT7 


* <P5  - P10>  ^ + <P5  - P9>  ^ + (P5  - V + <P>  - Pl5> 


5-15 


(3.9) 
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where  R . denotes  the  resistance  of  each  leakage  path  and  P denotes 
i-J  k 

the  pressure  of  the  leakage  source  or  sink.  The  leakage  flow  is  computed 
for  discrete  rotor  positions.  A curve  is  then  fit  to  the  computed  flow 
rates  to  obtain  the  leakage  as  a continuous  function  of  rotor  position. 

The  pressures  in  Equation  (3.9)  were  obtained  as  a function  of 
rotor  position  from  the  thermodynamic  model.  However,  the  pressure  of  a 
control  volume  is  dependent  on  the  leakage  flow  rate.  Hence,  it  was 
necessary  to  iterate  to  obtain  an  accurate  measure  of  the  flow.  This  was 
accomplished  by  initially  inserting  the  control  volume  pressures,  computed 
in  the  absence  of  leakage,  into  the  leakage  expressions.  The  leakage 
flow  rate  was  computed  and  then  input  into  the  thermodynamic  model  to 
recompute  the  control  volume  pressures.  Based  on  these  new  pressures, 
the  leakage  flow  was  recomputed.  This  process  was  repeated  until  the 
control  volume  pressures  varied  only  slightly  from  iteration  to  iteration. 


3.4  Heat  Transfer  Model 

3.4.1  One-Dimensional  Model.  A one-dimensional  analysis  of  the 
heat  transfer  between  the  working  fluid  and  the  stator  as  well  as  between 
the  working  fluid  and  the  end  plates  was  performed  to  obtain  as  estima- 
tion of  the  effect  of  heat  transfer  on  the  expander  performance.  The 
fluid  motion  over  the  stator  and  end  plates  was  modeled  as  flow  over  a 
flat  plate,  and  standard  flat  plate  relationships  were  used  to  compute 
the  convective  heat  transfer  coefficients.  The  electrical  analogy  to 
heat  transfer  was  employed  and  the  heat  transfer  rate  was  obtained  as  a 
function  of  the  various  thermal  resistances  (component  material,  boundary 
layer,  etc.)  and  the  temperature  difference  between  the  working  fluid  and 
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the  surroundings.  The  working  fluid  temperature  was  obtained  directly 
from  the  thermodynamic  model. 

The  simplified  heat  transfer  analysis,  while  permitting  a gross 
determination  of  the  heat  transfer  rates  and  their  effect  on  the  expander 
performance  relative  to  leakage  and  friction,  does  not  permit  an  accurate 
determination  of  the  temperature  profiles  in  the  stator,  rotor  or  end 
plates  (necessary  for  computations  of  thermal  expansion).  Furthermore, 
the  assumption  of  one-dimensional  heat  transfer  allows  only  approximate 
values  of  the  heat  transfer  rates  to  be  obtained. 

To  obtain  a better  representation  of  the  heat  transfer  in  the 
expander,  more  sophisticated  models  of  the  rotor,  stator  and  end  plates 
were  required.  The  approach  taken  in  modeling  these  components  will  be 
discussed  separately. 


3.4.2  Rotor,  Three-Dimensional  Model.  A finite  difference 
approach  was  used  to  approximate  the  heat  conduction  equation 


1 3T  = 9_T  1 3T  l 
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(3.10) 


in  the  expander  rotor.  Appendix  C illustrates  the  method  used  in  forming 
the  finite  difference  equations. 

Since  the  flow  pattern  in  the  control  volume  is  not  known,  it  is 
difficult  to  determine  the  heat  transfer  coefficient  between  the  steam 
and  the  rotor.  Therefore,  an  experimentally  measured  surface  temperature 
was  specified  as  the  boundary  condition  on  the  rotor  surface.  It  was 
assumed  that  for  a given  control  volume,  the  rotor  surface  temperature 


did  not  vary  in  the  axial  or  tangential  direction.  This  was  an 
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approximation  since  frictional  heat  generation  between  the  vane  and  rotor 
would  cause  a higher  surface  temperature  to  be  observed  near  the  vane 
slots. 

Due  to  friction,  heat  is  generated  at  the  vane  rotor  interface. 
Therefore,  the  boundary  condition  along  the  slot  walls  at  points  in 
contact  with  a vane  was  a specified  heat  flux.  The  total  amount  of  heat 
generated  was  determined  by  the  friction  model  described  earlier.  A 
fraction  of  the  heat  generated  is  transferred  into  the  rotor  and  the 
remainder  is  transferred  into  the  vane.  To  obtain  an  approximation  of 
the  fractions  of  the  heat  generated  going  into  the  vane  and  rotor,  an 
exact  solution  to  the  temperature  profiles  in  two  semi-infinite  solids 
in  frictional  contact  was  obtained.  The  temperatures,  at  the  surface  of 
contact,  for  the  solids  were  equated  resulting  in  the  following  expres- 
sions for  the  amount  of  heat  transferred  into  each  solid  (see  Appendix 
D)  : 


An  adiabatic  boundary  condition  was  assumed  for  that  portion  of 
the  vane  slot  not  in  contact  with  a vane. 

The  rotor  is  symmetric  in  the  axial  direction.  Therefore,  only 
one-half  of  the  axial  length  of  the  rotor  was  used  in  the  analysis.  The 


38 


3 


The  rotor  model  included  provisions  for  using  a specified  heat 
flux,  resulting  from  friction,  as  the  boundary  condition  for  the  rotor 
end  facing  the  end  plate.  However,  based  on  the  freedom  of  rotor  motion 
in  the  experimental  expander,  it  was  assumed  that  there  was  no  rotor-end 
plate  contact  and  hence,  no  resulting  frictional  energy  dissipation.  It 
was  also  assumed  that  the  viscous  energy  dissipation  between  the  oil  and 
rotor  only  resulted  in  the  increase  in  the  temperature  of  the  oil. 

Lack  of  end  plate-rotor  contact  eliminated  a major  source  of  heat 
generation.  With  this  source  of  heat  generation  eliminated,  the  ro>tor 
end  effects  were  neglected  and  the  heat  conduction  equation  solved  for 
two  dimensions.  Analysis  of  the  computed  rotor  temperature  profile 
revealed  that  due  to  the  rapid  small  amplitude  fluctuations  of  the  rot  - 
surface  temperature,  the  changing  surface  conditions  did  not  greatly 
effect  the  rotor  temperature  at  points  well  below  the  surface.  There- 
fore, to  permit  the  use  of  several  nodes  near  the  surface  without 
excessively  increasing  the  computation  time,  modifications  to  the 
analysis  were  made. 

The  mean  rotor  surface  temperature  was  computed  and  that  tempera- 
ture used  as  a boundary  condition  at  a depth  of  0.050  inches  below  the 
surface.  Utilizing  this  boundary  condition  permitted  considering  only 
that  portion  of  the  rotor  contained  between  two  vane  slots.  The 
boundary  condition  at  the  vane  slot  and  the  rotor  ends  were  not  changed. 
The  heat  transfer  to  the  rotor  was  computed  using  the  temperature 
gradient  near  the  surface.  Since  the  rotor  surface  temperature  did  not 
vary  greatly,  the  mean  rotor  surface  temperature  was  sufficiently  good 
for  use  in  computing  the  rotor  thermal  expansion. 
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3.4.3  Stator,  Three-Dimensional  Model.  To  obtain  the  temperature 
profiles  and  heat  transfer  rates,  a three-dimensional  analysis  of  the 
stator  was  required  due  to  the  location  of  the  inlet  and  exhaust  ports 
at  the  axial  center  of  the  stator.  The  ports  act  as  heat  sources  or 
sinks  causing  temperature  gradients  to  exist  between  the  ports  and  the 
ends  of  the  stator.  The  three-dimensional  steady  state  heat  conduction 
equation 
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was  transformed  into  finite  difference  form  for  the  stator  geometry,  as 
shown  in  Figures  10.  The  method  of  forming  the  finite 
difference  expressions  is  illustrated  in  Appendix  C. 

At  the  inside  surface,  heat  is  transferred  to  the  stator  by 
convection.  Additionally,  heat  is  generated  at  the  inside  surface  due 
to  friction  between  the  vanes  and  the  stator.  Therefore,  the  inside 
surface  boundary  condition  was 


h(T  - 
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gen 


(3.13) 


The  fluid  temperature,  T^  and  heat  generation  rate  seen  by  a 
particular  position  on  the  stator  surface  are  actually  functions  of  time. 
However,  to  permit  a steady  state  analysis,  the  time  averaged  wall 
temperature  and  heat  generation  rates  were  employed.  For  a given 


position  on  the  stator  surface,  the  surface  temperature  was  computed  as: 
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and  the  heat  generation  rate  was  computed  as: 


(3.14) 


(3.15) 


The  time  dependent  fluid  temperatures  and  frictional  heat 
generation  rates,  used  in  Equations  (3.14)  and  (3.15)  were  obtained  from 
the  thermodynamic  and  friction  models,  respectively. 

No  iteration  between  the  steam  temperatures  and  heat  transfer  rates 
was  done  as  in  the  case  of  the  steam  pressures  and  the  leakage  flow  rates 
since  the  heat  transfer  did  not  appear  to  affect  the  component  temperature 
profiles,  as  a result  of  changing  the  steam  temperature,  to  the  same 
degree  as  leakage  affected  the  control  volume  pressure.  The  amount  of 
frictional  heat  entering  the  stator  and  the  vane  was  computed  using 
Equation  (3.11). 

An  estimate  of  the  convective  heat  transfer  coefficient  was 
obtained  by  modeling  the  flow  over  the  stator  as  flow  over  a constant 
temperature  semi- inf inite  flat  plate.  The  mean  heat  transfer  coeffi- 
cients from  Reference  10  were  for 


LAMINAR  FLOW: 
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MIXED  LAMINAR  - TURBULENT  FLOW: 

h = £ Pr1/3(O.OJ6  Re  °‘8  - 836) 
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(3.17) 


Transition  from  laminar  to  turbulent  flow  was  based  on  a critical 
Reynolds  Number  of  5 x 10^ . The  characteristic  length  used  in  computing 
the  Reynolds  Number  and  heat  transfer  coefficients  was  choosen  as 


Lg  = (Ri  + HR)  ARC  . 


(3.18) 


The  velocity  used  in  computing  the  Reynolds  Number  was  computed  as 


vmean's  = (Ri  + HR)  W • 


(3.19) 


Since  the  outside  surface  of  the  stator  was  insulated,  an 
adiabatic  boundary  condition  was  employed  there. 

It  was  assumed  that  the  surfaces  of  the  inlet  and  exhaust  ports 
in  contact  with  the  working  fluid  were  at  the  inlet  and  exhaust 
temperature,  respectively. 

The  stator  was  symmetric  in  the  axial  direction  about  a plane 
cutting  through  the  center  of  the  inlet  and  exhaust  ports.  Therefore, 
only  one  half  of  the  axial  length  of  the  stator  was  used  in  computing 
the  temperature  profiles  and  heat  transfer  rates.  An  adiabatic  boundary 
condition  was  employed  at  the  axis  of  symmetry. 

In  order  to  avoid  a computer  storage  problem  which  would  have 
resulted  when  considering  the  stator  and  end  plates  together,  they  were 
analyzed  separately.  To  permit  separate  analysis,  it  was  assumed  that 


the  boundary  between  the  stator  and  end  plate  was  adiabatic. 
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3. A. 4 End  Plate,  Three-Dimensional  Model.  To  obtain  the 
temperature  profiles  and  heat  transfer  rates  for  the  end  plate.  Equation 
(3.12)  was  transformed  into  finite  difference  form.  The  transformation 
is  illustrated  in  Appendix  C.  The  end  plate  was  modeled  as  shown  in 
Figures  11. 

A combined  convection  and  heat  generation  boundary  condition  was 

employed  at  the  end  plate  surfaces  that  contact  the  working  fluid 

(-k  = h(T  - T ) + q"  ) . 

dZ  00  w gen 

Since  the  fluid  temperature  and  heat  generation  rate  varied  with 
time.  Equations  (3.1A)  and  (3.15)  were  used  to  compute  the  time  averaged 
fluid  temperature  and  heat  generation  rates  for  a particular  point  on 
the  end  plate.  The  time  dependent  fluid  temperatures  and  frictional 
heat  generation  rates  used  in  Equations  (3.1A)  and  (3.15)  were  obtained 
from  the  thermodynamic  and  friction  models,  respectively.  Like  the 
stator,  no  iteration  between  the  steam  temperature  and  heat  transfer 
rates  was  done.  The  fraction  of  the  total  frictional  heat  generated 
between  the  vane  and  end  plate  that  entered  the  end  plate  was  computed 
using  Equation  (3.11).  The  heat  transfer  coefficients  for  the  end 
plate  were  computed  using  Equations  (3.16)  and  (3.17).  The  character- 
istic length  used  in  computation  of  the  Reynolds  Numbers  and  heat 


transfer  coefficients  was  chosen  as: 
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The  velocity  used  in  computing  the  Reynolds  Number  was  computed  as: 
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Consistent  with  the  stator  analysis,  the  surface  of  the  end  plate 
which  contacted  the  stator  was  assumed  to  be  adiabatic  (—  = 0)  . 

t)Z 

Since  there  was  no  contact  between  the  rotor  and  end  plate,  there 
was  no  heat  generated  on  the  end  plate  surface  by  rotor  friction.  How- 
ever, there  was  leakage  between  the  rotor  and  end  plate.  To  model  this 
condition,  a convective  boundary  condition  was  employed  on  those 

surfaces  of  the  end  plate  which  were  covered  by  the  rotor  but  did  not 

Vl 

contact  a vane  (-  -r—  * — (T  - T ) ] . 
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The  periphery  of  the  hole  in  the  end  plate  through  which  the 

9T 

rotor  shaft  extended  was  also  assumed  to  be  an  adiabatic  surface  (~—  = 0)  . 

9r 

The  shaft  did  not  contact  the  end  plates  at  this  point,  so  no  frictional 
heat  was  generated.  Furthermore,  since  the  rotor  shaft  was  sealed  at 
both  ends,  it  was  assumed  that  steam  leakage  between  the  end  plate  hole 
and  rotor  shaft  would  be  small.  It  was  possible  that  natural  convection 
might  have  existed,  but  it  would  have  been  small  in  comparison  to  the 
other  modes  of  heat  transfer  to  the  end  plate. 

Since  the  end  plates  were  insulated,  the  outside  surfaces  of  the 
end  plates  were  assumed  adiabatic. 


3.5  Numerical  Procedure 

A finite  difference  method  was  used  to  obtain  solutions  to  the 
expander  heat  transfer  analysis. 
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3.5.1  Transient  Analysis.  The  rapidly  fluctuating  surface 
temperature,  used  as  a rotor  boundary  condition,  required  the  use  of 
small  time  steps  to  analytically  synthesize  the  temperature  fluctuations. 
The  stability  requirements  for  the  explicit  method  were  determined  using 
the  technique  of  Karplus  [11,12].  It  was  found  that  the  time  steps 
required  for  stability  were  not  prohibitive  in  view  of  the  rapidly 
changing  rotor  surface  temperature.  Therefore,  the  explicit  method  was 
utilized  in  the  rotor  heat  transfer  analysis. 

2 * 

When  employing  the  explicit  method,  the  Laplacian  (V  T)  was  put 

into  finite  difference  form  using  the  method  of  Taylor's  Series 

expansion  [11,13].  A forward  difference  approximation  was  used  to 

9T 

approximate  the  transient  term  (-^-) . This  is  illustrated  in  Appendix  C. 

3.5.2  Steady  State  Analysis.  The  method  of  Taylor's  Series 

expansion  was  used  to  obtain  the  finite  difference  approximation  to  the 
2 

Laplacian  (V  T).  The  set  of  algebraic  equations  that  resulted  from  the 
finite  difference  analysis  were  solved  using  the  iterative  method  of 
Gauss-Seidel  with  the  incorporation  of  Young's  over  relaxation  factor 
parameter  [7]. 
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CHAPTER  IV 

EXPERIMENTAL  AND  THEORETICAL  RESULTS 

4 . 1 Experimental  Conditions 

The  expander,  described  in  Chapter  2,  was  tested  with  a 22.5°  and 

a 45°  arc  of  admission.  For  each  arc  of  admission,  supply  pressures  of 

150  psia  and  115  psia  were  employed.  The  supply  temperatures  for  the 

150  psia  and  115  psia  supply  pressures  were  370°F  + 10°F  and  360°F  + 10°F, 

respectively.  The  speed  was  varied  from  840  to  1913  RPM. 

Initial  testing  of  the  expander  was  conducted  using  no  internal 

lubrication.  However,  in  an  attempt  to  reduce  component  wear  and  to 

decrease  leakage  by  partially  blocking  the  leakage  paths  with  oil,  all 

the  data  included  herein  was  collected  using  an  oil  flow  rate  of  0.05 

lb  /min.  This  compared  with  a steam  flow  rate  on  the  order  of  5 lb  /min. 
m m 

4.2  Experimental  Results 

4.2.1  Power . Figures  12  and  13  are  representative  pressure 
versus  volume  diagrams  constructed  from  the  measured  internal  expander 
pressures.  The  experimentally  determined  data  points  are  indicated  on 
the  curves.  However,  portions  of  the  diagram  were  constructed  by 
assuming  the  pressures  in  the  control  volume  at  the  inlet  and  the  exhaust 
ports  were  the  supply  and  exhaust  pressures,  respectively.  This  is  a 
good  assumption  since  the  vane  expander  is  an  easy  breathing  device. 

This  has  been  verified  analytically  in  Reference  1,  and  by  the 


P* (DIMENSIONLESS  PRESSURE) 


48 


0.0  0.05  0.10  0.15  0.20  0.25  0.30 

V*  (DIMENSIONLESS  VOLUME) 


Figure  12. 


Pressure  Versus  Volume  for  150  psia  Supply 
Pressure  22.5°  Arc  of  Admission 
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V*  (DIMENSIONLESS  VOLUME) 


Figure  13.  Pressure  Versus  Volume  for  150  psia  Supply 
Pressure  45°  Arc  of  Admission 
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experimental  work  discussed  in  Reference  4.  Included  in  Figures  12  and 
13  are  the  ideal  P*  -V*  curves  for  the  expander.  The  ideal  curve  is 
based  on  the  analysis  of  an  expander  with  no  friction,  leakage  or  heat 
transfer.  During  the  expansion  portion  of  the  cycle,  the  measured  control 
volume  pressures  are  higher  than  the  ideal  values.  This  is  primarilly 
due  to  leakage  from  the  region  exposed  to  the  inlet  port  into  the  lower 
pressure  control  volumes.  The  rapid  decrease  in  control  volume  pressure 
prior  to  opening  the  exhaust  port  is  due  to  leakage  out  of  the  expanding 
control  volume  to  the  low  pressure  region  exposed  to  the  exhaust  port. 

The  increased  area  of  the  actual  P*  -V*  diagram,  compared  to  the 
ideal  diagram,  illustrates  that  leakage  increased  the  expander  power 
output . 

The  indicated  work  per  revolution  done  by  the  expander  was 
numerically  computed  as 


W 


PdV 

rev 


(4.1) 


using  the  measured  pressures  and  the  known  volume  variation.  Figures  14 
and  15  show  the  indicated  power  as  a function  of  speed  for  the  22.5°  and 
45° arcs  of  admission  and  supply  pressures  of  150  psia  and  1J5  psia.  Due 
to  the  larger  arc  of  admission,  higher  indicated  power  outputs  were 
obtained  for  the  45°  arc  of  admission.  The  difference  in  indicated 
power,  for  the  two  arcs  of  admission,  would  have  been  greater  had  the 
back  pressure  on  the  expander  not  increased,  due  to  increased  mass  flow, 
at  the  larger  arc  of  admission.  The  shaft  power,  as  measured  by  the 
dynamometer,  and  the  frictional  power  loss  are  also  shown  in  Figures  14 


Figure  14.  Power  Versus  Speed  for  150  psia  Supply  Pressure 
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Figure  15.  Power  Versus  Speed  for  115  psia  Supply  Pressure 
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and  15.  The  frictional  power  loss  was  determined  as  the  difference 
between  the  indicated  power  and  the  measured  shaft  power.  Over  the 
entire  speed  range  for  150  psia  supply  pressure  and  over  the  majority 
of  the  speed  range  for  115  psia  supply  pressure,  the  45°  arc  of 
admission  data  shows  a greater  frictional  power  loss.  This  was  due  to 
higher  pressure  forces  being  applied  to  the  vane  over  a greater  portion 
of  the  cycle  for  the  larger  arc  of  admission.  The  larger  pressure 
forces  produce  an  increase  in  vane  friction.  Figure  15  shows  that  for 
a supply  pressure  of  115  psia,  the  frictional  power  loss  for  the  22.5° 
arc  of  admission  exceeded  the  frictional  power  loss  for  the  45°  arc  of 
admission  at  speeds  above  1650  RPM.  This  was  probably  due  to  a 
variation  in  the  coefficient  of  sliding  friction.  However,  since  the 
coefficient  of  friction  is  dependent  on  many  factors  such  as  load,  and 
the  amount  and  location  of  lubricating  oil  and  liquid  water,  it  is 
difficult  to  explain  why  the  variation  may  have  occurred. 

Figure  14  shows  that  although  the  indicated  and  frictional  power 
increased  linearly  with  speed,  the  shaft  power  curve  begins  to  bend  near 
1500  RPM,  for  a supply  pressure  of  150  psia.  This  is  due  to  a decrease 


in  the  mechanical  efficiency  (mechanical  efficiency  = 


shaft  horsepower  . 
indicated  horsepower"1 


at  the  higher  speed.  The  decrease  in  mechanical  efficiency  is  a result 
of  the  frictional  power  loss  becoming  a larger  percentage  of  the  indicated 
power.  The  variation  of  mechanical  efficiency  with  speed  for  a supply 
pressure  of  150  psia  is  shown  in  Figure  16. 

For  a supply  pressure  of  115  psia,  and  a 22.5°  arc  of  admission, 
the  shaft  power  curve  begins  to  bend  at  approximately  1300  RPM.  For  this 


supply  pressure,  the  45°  arc  of  admission  data  show  a continuous  increase 
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in  shaft  power  with  speed.  Figure  17  shows  that  the  mechanical 
efficiency,  for  a 22.5°  arc  of  admission  and  115  psia  supply  pressure, 
drops  off  quite  rapidly.  This  decrease  in  mechanical  efficiency,  caused 
by  friction,  results  in  the  bending  of  the  22.5°  arc  of  admission  shaft- 
power  curve.  The  mechanical  efficiency  for  the  45°  arc  of  admission  and 
115  psia  supply  data,  while  varying  with  speed,  shows  no  definite  trend 
downward  at  the  higher  speeds. 

4.2.2  Mass  Flow  Rate.  Figures  18,  19,  20  and  21  show  the 
measured  mass  flow  rate,  the  ideal  flow  rate  and  the  leakage  flow  rate 
as  a function  of  speed  for  45°  and  22.5“  arcs  of  admission  and  115  psia 
and  150  psia  supply  pressures.  In  all  cases,  the  flow  rate  of  steam 
increased  with  speed.  Ideally,  for  an  easy  breathing  expander,  the 
control  volume  contains  the  same  mass  of  steam  per  revolution,  independ- 
ent of  speed.  Therefore,  as  the  speed  increases  the  mass  flow  rate 
increases.  The  ideal  mass  flow  rate  was  computed  as  the  product  of  the 
number  of  control  volumes,  the  volume  of  a control  volume  at  intake 
closing,  the  density  of  the  supply  steam,  and  the  speed  in  revolutions 
per  minute.  The  leakage  flow  was  computed  as  the  difference  between  the 
measured  flow  rate  and  the  ideal  flow  rate. 

In  all  cases  studied,  the  leakage  flow  remained  approximately 
constant  with  speed.  The  clearances  in  the  expander  do  not  change  with 
speed  and  for  an  easy  breathing  expander,  the  pressures  in  the  control 
volume  do  not  vary  greatly  with  speed.  Since  the  leakage  driving 
potential  and  flow  resistance  do  not  vary  with  speed,  it  is  consistent 
that  the  leakage  flow  rate  would  not  change  greatly  with  speed. 


MASS  FLOW  RATE  (lb„/min) 


Figure  19.  Mass  Flow  Rate  Data  for  150  psia  Supply  Pressure 
45°  Arc  of  Admission 


MASS  FLOW  RATE  (lb„/min.) 


Figure  21.  Mass  Flow  Rate  Data  for  115  psia  Supply  Pressure 
45°  Arc  of  Admission 
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Comparison  of  Figures  18  and  19  indicates  that  for  the  same 
supply  pressure,  the  leakage  flow  increases  with  arc  of  admission.  The 
same  trend  is  observed  by  comparison  of  Figures  20  and  21.  The 
resistance  to  leakage  flow,  for  a control  volume  at  any  specified 
position,  does  not  change  with  arc  of  admission;  but  the  pressures  in 
the  control  volumes  are  maintained  at  higher  values  over  a greater 
portion  of  the  cycle  with  the  larger  arc  of  admission.  This  greater 
driving  potential  is  apparently  the  cause  for  higher  leakage  at  the 
larger  arc  of  admission. 

Comparisons  of  Figures  18  and  20,  and  Figures  19  and  21,  show 
that  for  the  same  arc  of  admission,  the  leakage  flow  increases  with 
supply  pressure.  This  would  also  appear  to  be  due  to  the  increased 
driving  potential. 

4.2.3  Adiabatic  Expansion  Efficiency.  Figure  22  shows  the 
adiabatic  expansion  efficiency  for  both  arcs  of  admission  at  a supply 
pressure  of  150  psia.  The  adiabatic  expansion  efficiency  is  defined  as 
the  ratio  of  the  work  per  pound  of  steam,  to  the  isentropic  enthalpy 
change  from  the  supply  conditions  to  the  exhaust  manifold  pressure.  The 
low  values  of  efficiency,  shown  in  Figure  22,  are  primarilly  the  result 
of  leakage  and  friction. 

Leakage  has  its  greatest  effect  at  low  speeds.  If  all  losses  but 
leakage  are  neglected,  and  it  Is  assumed  that  leakage  does  not  increase 
the  power  output,  the  adiabatic  expansion  efficiency  can  be  written  as 

^ADIABATIC  1 ~ • (A‘2) 

EXPANSION  mT 
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where  is  the  leakage  flow  rate  and  m^  is  the  total  flow  rate 

leakage  flow  increases,  the  efficiency  decreases.  The  effect  of  leakage 
on  the  adiabatic  expansion  efficiency  is  reduced  somewhat  since  leakage 
does  increase  the  expander  power  output.  Figures  18  through  21  showed 
that  the  experimentally  computed  leakage  flow  rates  did  not  vary  greatly 
with  speed.  However,  the  ideal  flow  rate  is  directly  proportional  to 
speed.  This  implies  that  as  the  speed  increases,  the  ratio  of  leakage 
flow  to  total  flow  decreases  resulting  in  a smaller  adiabatic  expansion 
efficiency  reduction  due  to  leakage. 

The  effect  of  friction  on  the  adiabatic  expansion  efficiency 
increases  with  increased  speed.  For  an  easy  breathing  expander,  the 
frictional  forces  increase  slightly  with  speed  due  to  larger  vane 
accelerations.  The  distance  per  unit  time  over  which  the  frictional 
forces  act  also  increases  directly  with  speed.  Therefore,  as  shown  in 
Figures  14  and  15,  the  frictional  power  loss  increases  with  speed. 

Under  the  influence  of  friction  and  leakage,  a peak  in  the 
adiabatic  expansion  efficiency  versus  speed  will  occur.  Low 
efficiencies  at  speeds  less  than  the  peak  efficiency  speed  would  be 
primarilly  due  to  leakage.  The  low  efficiencies  at  speeds  greater  than 
peak  efficiency  speed  would  be  primarilly  due  to  friction. 

Figure  22  shows  that  at  15U  psia  supply  pre  ure  and  arcs  of 
admission  of  22.5°  and  45°,  the  experimentally  determined  efficiency 
continuously  increased  with  speed.  This  was  due  to  a decrease  in  the 
ratio  of  leakage  flow  rate  to  total  mass  flow  rate.  For  the  conditions 
shown  in  Figure  22,  the  peak  adiabatic  expansion  efficiency  was  not 
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Figure  23  shows  the  same  trend  for  the  45°  arc  of  admission  and 
115  psia  supply  pressure.  However,  for  the  22.5°  arc  of  admission  and 
115  psia  supply  pressure,  a peak  in  the  experimental  adiabatic  expansion 
efficiency  was  observed  at  approximately  1450  RPM.  The  decrease  in 
efficiency  with  speed  above  1450  RPM  was  due  to  frictional  power  losses. 
The  severe  effect  of  friction  for  a supply  pressure  of  115  psia  and  an 
arc  of  admission  of  22.5°  was  previously  illustrated  by  the  rapid  drop 
in  mechanical  efficiency  shown  in  Figure  17. 

4.2.4  Temperature  Measurement  Results.  The  liner  and  end  plate 
temperatures  were  measured  for  determining  thermal  expansions  and  for 
comparison  with  the  temperatures  predicted  by  the  heat  transfer  analysis. 
The  rotor  surface  temperature  was  measured  so  that  it  could  be  used  as  a 
boundary  condition  in  the  heat  transfer  analysis.  It  also  permitted 
computation  of  the  approximate  rotor  thermal  expansion. 

Figures  24  through  27  show  the  measured  stator  temperature 
profiles  (the  points  indicate  the  experimental  data)  for  the  22.5°  and 
45°  arcs  of  admission  and  supply  pressures  of  150  psia  and  115  psia. 

The  temperatures  were  measured  0.0625  inches  from  the  expander  end  plate 
and  0 . 050  inches  f rom  the  inside  surface  of  the  stator.  The  temperatures 
have  been  normalized  by  the  supply  temperature.  Although  the  stator 
temperature  is  higher  on  the  inlet  side  of  the  expander,  the  temperature 
profile  is  not  severely  distorted.  Temperature  differences  between  the 
inlet  and  exhaust  side  of  the  expander  can  result  in  distortion  of  the 
expander  stator.  Distortion  of  the  stator  would  make  it  more  difficult 
to  maintain  small  clearances  between  the  rotor  and  end  plates  without 
producing  large  frictional  losses. 
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Comparison  of  Figures  24  through  27  shows  that  there  is  no  i 

{ 

J 

, significant  difference  between  the  stator  temperature  profiles  for  the 

22.5°  and  the  45°  arc  of  admission  cases.  Also,  the  supply  pressure 
does  not  significantly  affect  the  normalized  temperature  profile. 

Figures  28  through  31  illustrate  that  there  is  no  severe 
I distortion  of  the  end  plate  temperature  profile.  The  arc  of  admission 

and  supply  pressure  also  does  not  greatly  affect  the  end  plate 
normalized  temperature  profile. 

Figure  32  shows  representative  curves  of  the  normalized  rotor 
surface  temperature  as  a function  of  rotor  position.  The  data 
indicates  that  the  rotor  surface  temperature  increases  when  the  surface 
is  exposed  to  the  high  temperature  supply  steam.  The  surface  tempera- 

i ! 

ture  begins  to  decrease  when  the  inlet  port  is  closed  and  continues  to 
decrease  until  the  inlet  reopens. 

The  variation  of  the  rotor  surface  temperature  with  rotor  position 
\ is  not  large.  This  is  due  to  the  small  time  during  each  revolution  for 

the  transfer  of  heat  between  the  rotor  and  the  working  fluid. 

The  rotor  surface  temperature  at  the  point  of  the  inlet  opening 
was  below  the  saturation  temperature  of  the  supply  steam.  This  would 
indicate  that  some  condensation  of  the  steam  occurred.  Condensation 
will  be  discussed  in  a later  section. 

4.2.5  Expander  Heat  Loss.  The  net  heat  loss  from  the  expander 
was  determined  by  applying  the  first  law  of  thermodynamics  to  the 
expander.  Table  2 summarizes  the  computed  heat  losses  for  the  tests 
conducted.  The  values  predominantly  range  from  3000  to  7500  BTU/hr 
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45°  ARC  OF  ADMISSION 
150  PSIA  SUPPLY  PRESSURE 


Figure  29.  End  Plate  Temperature  Profile  Data  for  150 
p8ia  Supply  Pressure  45°  Arc  of  Admission 
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ROTOR  POSITION  (RADIANS) 


Figure  32.  Rotor  Surface  Temperature  Data 
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TABLE  2 
Expander  Heat 

Loss  Data 

Supply  Pressure 
(psia) 

Arc  of  Admission 
(Degrees) 

Speed  (RPM) 

Heat  Loss  <BTU) 
hr 

150 

22.5 

1007 

3639 

150 

22.5 

1194 

605 

150 

22.5 

1488 

3219 

115 

22.5 

1083 

4942 

115 

22.5 

1109 

3186 

115 

22.5 

1350 

4024 

115 

22.5 

1693 

4486 

150 

45 

840 

6661 

150 

45 

1138 

6619 

150 

45 

1369 

7188 

150 

45 

1479 

4581 

150 

45 

1539 

6072 

150 

• 45 

1906 

6359 

115 

45 

895 

7481 

115 

45 

1285 

4081 

115 

45 

1296 

2801 

115 

45 

1477 

6028 
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(1.18  to  2.95  horsepower).  This  indicates  that  the  heat  loss  was  of  the 
same  order  of  magnitude  as  the  power  output.  The  variations  observed  in 
the  expander  heat  loss  are  believed  to  be  primarilly  due  to  the  error 
sustained  in  measuring  the  inlet  and  exhaust  enthalpy.  An  error  of  25 
BTU/lb^  in  determining  the  inlet  and  exhaust  enthalpy  difference  results 
in  an  error  in  the  heat  loss  of  the  order  of  4500  BTU/hr.  This  enthalpy 
difference  of  25  BTU/lbm  is  approximately  27.  of  the  inlet  enthalpy. 

Estimates  of  the  magnitude  of  the  heat  loss  to  the  slip  ring  and 
pressure  transducer  cooling  adapters  were  of  the  same  order  as  the 
experimentally  determined  heat  losses.  This  source  of  energy  loss  would 
be  eliminated  in  an  operational  expander. 

4 . 3 Comparisons  with  Theoretical  Model 

4.3.1  Mass  Flow  Rate.  Computation  of  the  mass  flow  rate  depends 
on  the  determination  of  the  internal  expander  clearances.  The  rotor-end 
plate  and  vane-end  plate  clearances  were  computed  based  on  the  initial 
cold  clearances  and  the  relative  thermal  expansion  of  the  components. 

The  cold  clearances  were  determined  from  micrometer  measurements 
taken  during  assembly  of  the  expander.  The  uncertainty  in  the  determin- 
ation of  the  cold  clearances  was  estimated  to  be  0.001  inches. 

The  thermal  expansions  were  computed  using  the  measured  component 
temperatures.  Although  the  thermal  expansions  of  the  components  varied 
with  the  supply  and  exhaust  temperatures,  the  variations  were  small 
relative  to  the  accuracy  of  the  cold  clearance  measurements.  The  rotor- 
end  plate  clearance,  under  operating  conditions,  was  calculated  to  be 
0.003  inches.  The  vane-end  plate  clearance  was  calculated  to  be  0,005 
inches.  The  larger  clearance  between  the  vanes  and  the  end  plates  was 
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due  to  the  small  coefficient  of  the  thermal  expansion  of  the  vanes 
2.7  x 10  ^ in/in°F,  compared  to  5.6  x 10  ^ in/in°F  for  the  rotor. 

The  computed  rotor-end  plate  and  vane-end  plate  clearances  were 
used  in  the  leakage  model  to  compute  the  leakage  flow  rate.  The  vane-end 
plate  clearance  was  adjusted  from  0.005  inches  to  0.006  inches  due  to 
wear  of  the  carbon  vanes  and  end  plates.  It  also  permitted  better 
agreement  to  be  obtained  with  the  experimentally  determined  leakage  flow 
rates.  Adjustment  of  the  clearance  by  0.001  inches  was  felt  to  be 
justifiable  since  the  change  was  within  the  limits  of  the  cold  clearance 
measurements. 

In  computing  the  theoretical  leakage  flow  rates,  all  leakage  paths 
except  those  between  the  end  plate  and  the  rotor  and  the  end  plate  and 
the  vanes  were  neglected.  Figure  33a  shows  a comparison  between  the 
analytically  and  experimentally  determined  leakage  flow  rates  for  a 22.5° 
arc  of  admission  and  150  psia  supply  pressure.  The  predicted  leakage 
flow  rates  were  16%  below  the  measured  values.  It  is  believed  the  error 
is  due  to  the  assumption  that  there  was  no  leakage  between  the  vane  tip 
and  the  stator.  Clearances  between  the  vane  and  the  stator  may  have 
occurred  as  a result  of  the  vane  being  pushed  into  its  slot  by  high 
pressure  supply  steam.  Clearances  between  the  vane  and  stator  may  also 
have  occurred  due  to  vane  chatter  resulting  from  the  stator  being  rough, 
or  "out  of  round".  It  is  possible  that  any  vane-stator  clearance  would 
result  in  more  leakage  from  the  regions  exposed  to  supply  pressure  to 


those  control  volumes  in  which  the  closed  expansion  process  was  occurring. 


Since  vane-stator  clearances,  resulting  from  vane  chatter,  are  a 
function  of  time  and  rotor  position,  a reasonable  determination  of  the 
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Figure  33a.  Leakage  Flow  Rate  Versus  Speed 
Figure  33b.  Horsepower  Versus  Speed 
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vane  stator  clearance  was  not  possible  and  no  attempts  were  made  to  Include 
leakage  between  the  vane  and  stator. 

4.3.2  Power  Output.  The  analytically  determined  leakage  flow 
rates  were  input  into  the  thermodynamic  model  to  determine  the  effect  of 
leakage  on  the  expander  power  output. 

Figure  33b  shows  the  experimentally  determined  indicated  power 
output,  the  predicted  power  output  with  leakage  affects  included,  and 
the  ideal  power  output.  The  ideal  power  output  was  computed  assuming  no 
leakage,  heat  transfer  or  friction. 

The  analytical  power  output,  computed  considering  only  leakage 
effects,  is  7.7 % below  the  measured  indicated  power  at  1007  RPM.  This 
difference  Increases  to  30.4%  at  1820  RPM.  The  difference  between  the 
predicted  power  output  and  the  measured  indicated  power  output  is 
apparently  the  result  of  the  leakage  model  not  predicting  the  correct 
amount  of  leakage.  As  previously  discussed,  the  leakage  was  computed 
assuming  a perfect  seal  between  the  vane  and  the  stator.  Leakage  between 
the  vane  and  the  stator  would  result  in  more  leakage  from  those  regions 
at  supply  pressure  to  the  control  volumes  undergoing  the  closed  expansion 
process.  Leakage  of  this  nature  would  result  in  increased  expander  power 
output. 

The  effect  of  heat  transfer  on  power  was  analytically  computed  by 
using  the  theoretically  determined  heat  transfer  rates  in  the  thermo- 
dynamic model.  The  theoretical  indicated  power  was  computed  in  the 
presence  of  both  leakage  and  heat  transfer  for  a 22.5°  arc  of  admission 
and  ISO  psia  supply  pressure.  Heat  transfer  reduced  the  theoretical 
power  output  by  14%,  from  2.40  horsepower  to  2.06  horsepower,  at  1007  RPM. 
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At  1820  RPM,  heat  transfer  reduced  the  power  by  13.5%,  from  3.25  to 
2.81  horsepower. 

Comparison  between  the  analytical  (leakage  and  heat  transfer 
included)  and  experimental  indicated  power  outputs  for  a 22.5°  arc  of 
admission  and  150  psia  supply  pressure  shows  an  error  of  20.7%  at  1007 
RPM  and  39%  at  1820  RPM.  It  is  felt  that  the  error  is  predominately 
the  result  of  the  inability  to  accurately  determine  the  effects  of 
leakage  on  the  expander  power  output. 

Figure  34  shows  a comparison  between  the  analytical  and  experi- 
mental frictional  power  loss  for  a supply  pressure  of  150  psia  and  an 
arc  of  admission  of  22.5°. 

The  theoretical  and  experimental  curves  are  identical  within  the 
accuracy  of  the  plot.  Good  agreement  was  also  obtained  for  the  45°  arc 
of  admission  and  150  psia  supply  steam.  At  1479  RPM,  and  a 45°  arc  of 
admission,  the  experimental  frictional  power  loss  was  1.69  horsepower 
while  the  theoretical  value  was  1.71  horsepower.  At  1913  RPM,  and  a 
45°  arc  of  admission,  the  experimental  value  was  2.26  horsepower  and  the 
theoretical  value  2.25  horsepower. 

The  coefficient  of  friction  used  in  the  model,  for  all  surfaces, 
was  0.075.  This  value  was  obtained  by  inputting  the  measured  steam 
pressures  into  the  friction  model  and  adjusting  the  coefficient  of 
friction  until  the  predictions  matched  the  experimental  data.  The 
theoretical  frictional  power  loss  was  computed  using  the  coefficient  of 
friction  of  0.075  and  the  steam  pressures  computed  by  the  thermodynamic 
model,  with  leakage  effects  included.  A test  of  the  validity  of  the 


friction  model  was  that  a single  coefficient  of  friction  allowed  accurate 


prediction  of  the  frictional  power  loss  at  all  speeds  and  supply  condi- 
tions except  for  the  115  psia  supply  pressure  and  22.5s  arc  of  admission 
case.  However,  the  trends  in  the  power  output  and  efficiency  data  for 
this  latter  case  appear  to  indicate  an  actual  variation  in  the  coeffi- 
cient of  friction. 


4.3.3  Temperature  Profiles  and  Heat  Transfer  Coefficient.  A 
comparison  of  the  experimental  and  theoretical  stator  temperature 
profiles  is  given  in  Figure  35.  The  maximum  error  is  13%.  This  error 
was  considered  to  be  reasonably  good  due  to  the  complex  geometry  and 
lack  of  information  regarding  the  heat  transfer  coefficient  between  the 
steam  and  the  stator. 

The  experimental  and  theoretical  end  plate  temperature  profiles 
are  compared  in  Figure  36.  The  maximum  error  is  8%. 

The  measured  rotor  surface  temperature  was  used  as  a boundary 
condition  in  computing  the  heat  transfer  to  the  rotor.  The  computed 
rotor  heat  transfer  was  used  in  conjunction  with  the  measured  surface 
temperature  and  inlet  steam  temperature  to  estimate  the  heat  transfer 
coefficient  between  the  steam  and  the  rotor.  The  values  obtained  are 
only  approximate  since  at  any  instant  in  time  the  steam  temperature  is 
not  accurately  known.  However,  at  the  inlet  the  steam  was  assumed  to  be 
at  the  supply  temperature.  Under  this  assumption,  the  computed  heat 
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transfer  coefficient,  at  a speed  of  1007  RPM,  ranged  from  610  BTU/hrft  #F 
when  the  rotor  was  initially  exposed  to  the  supply  steam,  to  approximately 
14,000  BTU/hrft  SF,  Just  before  the  inlet  was  closed.  The  variation  in 
the  heat  transfer  coefficient  apparently  resulted  from  the  turbulence 


and  steam  impingement  on  the  rotor  surface.  The  measured  rotor  surface 
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temperatures  used  in  computing  the  rotor  heat  transfer  were  obtained  at 
the  circumferential  center  of  the  control  volumes.  When  the  rotor 
surface  of  a particular  control  volume  is  first  exposed  to  the  inlet 
port,  the  center  of  the  control  volume  is  restricted  from  seeing  the 
high  temperature  steam  by  the  small  rotor-stator  clearance  which  exists 
midway  between  the  inlet  and  exhaust  ports.  Since  the  center  of  the 
control  volume  would  not  be  exposed  to  direct  impingement  by  the  supply 
steam  and  the  turbulence  created  by  the  filling  process,  the  calculated 
heat  transfer  coefficient  would  be  expected  to  be  low  near  inlet  opening. 
However,  when  the  center  of  the  control  volume  becomes  exposed  to  the 
incoming  steam,  impingement  of  the  steam  on  the  center  of  the  rotor 
coupled  with  turbulence  created  by  the  filling  process,  would  cause  the 
heat  transfer  coefficient  to  be  quite  high. 

The  heat  transfer  coefficients  at  the  exhaust  ranged  from 
420  BTU/hr  ft^°F  when  the  exhaust  opened  to  a maximum  of  2800  BTU/hr  ft^°F 
as  the  exhaust  closing  was  approached.  The  increase  in  the  heat  transfer 
coefficient  during  the  exhaust  process  was  probably  due  to  an  increase  in 
the  turbulence  of  the  steam  as  it  was  forced  out  of  the  control  volume  by 
the  reduction  in  the  volume  of  the  control  volume. 

These  heat  transfer  coefficients  are  only  approximate  values; 
however  they  do  provide  a starting  point  for  computation  of  the  rotor 
heat  transfer  without  knowledge  of  the  rotor  surface  temperature.  It 
should  be  pointed  out  that  any  computation  of  the  instantaneous  heat 
transfer  in  the  expander  using  the  working  fluid  temperature  and  a heat 
transfer  coefficient  is  only  a rough  approximation.  This  is  due,  not 
only  to  the  uncertainty  of  the  heat  transfer  coefficient  but  also  to  an 
uncertainty  in  the  knowledge  of  the  instantaneous  steam  temperature. 
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Figure  37  illustrates  the  variation  in  the  ratio  of  the  rotor 
surface  temperature  to  the  steam  saturation  temperature  as  a function  of 
rotor  position.  The  data  indicate  that  during  a small  portion  of  each 
revolution,  primarilly  during  the  filling  process,  the  rotor  surface 
temperature  is  below  the  steam  saturation  temperature.  This  condition 
would  permit  condensation  of  some  steam  during  filling  of  the  control 
volume. 

An  order  of  magnitude  analysis  of  the  degree  of  condensation  can 
be  made  by  assuming  that  the  amount  of  steam  that  could  be  condensed 
during  the  filling  process  can  be  approximately  written  as 


h ^SlS^SAT  " TROTOR)At 


(4.1) 


where  h is  an  average  convective  heat  transfer  coefficient  and  h^  is 
the  latent  heat  of  vaporization  of  the  steam.  The  convective  heat 
transfer  coefficient  for  the  condensation  process  is 


k/6  , 


(4.2) 


where  k is  the  thermal  conductivity  of  liquid  water  and  6 is  the 
thickness  of  the  condensed  liquid  film. 

From  geometric  considerations,  the  amount  of  liquid  condensed 
during  the  filling  process  can  be  approximately  written  as 
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where  (R^A8Z)  is  the  rotor  surface  area  of  one  control  volume.  Equating 
Equations  (4.1)  and  (4.3)  and  using  Equation  (4.2)  gives 
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Solving  Equation  (4.4)  for  h and  substitution  of  h into  Equation  (4.1) 
permits  an  estimate  of  the  amount  of  steam  condensed  in  each  control 
volume  during  the  filling  process.  The  approximate  mass  flow  rate  of 
condensed  steam  is 
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Based  on  this  calculation,  condensation  rates  of  0.279  lb  /min  at 

m 

1820  RPM  and  0.207  lb  /min  at  1007  RPM  have  been  calculated. 

m 

The  heat  transfer  coefficients  required  to  condense  the  steam  are 
25,789  BTU/hr  ft2°F  at  1820  RPM  and  19,209  BTU/hr  ft2°F  at  1007  RPM. 


Lower  heat  transfer  coefficients  would  result  in  less  condensation,  at  a 
particular  speed,  due  to  insufficient  time  to  condense  the  steam.  The 
rotor  heat  transfer  coefficients  computed  from  the  rotor  heat  transfer 

analysis,  at  1007  RPM,  have  a minimum  value  at  the  inlet  of  610  BTU/hr 

2 2 
ft  °F  and  a maximum  of  14,000  BTU/hr  ft  °F.  Since  these  values  are  less 

than  that  computed  by  Equation  (4.4),  the  actual  amount  of  steam  con- 
densed on  the  rotor  would  be  less  than  the  calculated  values. 

It  is  believed  that  condensation  on  the  stator  and  the  end  plates 
would  be  less  than  on  the  rotor  since  the  surfaces  of  these  components 
do  not  see  the  relatively  cold  exhaust  steam  as  long  as  the  rotor  surface 
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does.  No  stator  or  end  plate  surface  temperatures  were  measured.  The 
measured  end  plate  and  stator  temperatures  at  0.065  inches  below  the 
surface  are  slightly  less  than  the  supply  temperature  at  the  inlet. 
However,  the  temperature  at  the  surface  of  these  components  would  be 
higher  due  to  convection  from  the  supply  steam  and  frictional  heat 
generation  at  the  surface. 

Since  the  predicted  condensation  rates  are  less  than  10%  of  the 
difference  between  the  measured  flow  rate  and  the  ideal  flow  rate,  it 
can  be  concluded  that  the  difference  was  primarilly  due  to  leakage.  It 
can  also  be  concluded  that  the  difference  between  the  predicted  power 
output  (computed  with  leakage  affects  only)  and  the  experimentally 
determined  indicated  po*rer  shown  in  Figure  33b  is  not  due  to  condensation 
and  re-evaporation.  The  estimated  condensation  rate  is  not  of  sufficient 
magnitude  to  cause  the  difference. 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 

5.1  Summary 

A commercially  available  air  expander  modified  to  operate  on 

steam  was  tested  using  150  psla  and  115  psia  supply  steam.  The  speed 

was  varied  from  840  to  1913  RPM.  Under  these  operating  conditions,  the 

power  output  ranged  from  0.5  horsepower  to  2.85  horsepower.  Mechanical 

efficiencies  ranged  from  65%  to  40%.  The  experimentally  determined 

adiabatic  expansion  efficiencies  varied  from  9.5%  at  the  low  speeds  to 

20%  at  the  high  speeds.  Leakage  flow  rates,  ranging  from  2.40  lb  /min 

to  3.90  lb  /min, were  obtained.  The  leakage  flow  increased  with 
m 

increased  arc  of  admission  and  supply  pressure.  The  total  mass  flow  of 

the  expander  varied  from  3.1  lb  /min  to  5.8  lb  /min. 

m m 

The  expander  was  analytically  modeled  to  permit  the  determination 
of  the  expander  performance  and  operating  component  temperatures.  j 

The  models  predicted  the  frictional  power  loss  very  well.  The 
leakage  flow  rates  were  predicted  with  an  error  of  16%.  The  heat 
transfer  models  predicted  the  stator  and  end  plate  temperature  profiles 
with  an  accuracy  on  the  order  of  10%.  Indicated  power  predictions  were 
computed  with  errors  ranging  from  20.7%  at  the  low  speeds  to  39%  at  the 
high  speeds. 
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5.2  Conclusions 

1)  The  adiabatic  expansion  efficiency  is  severely  reduced  by 
friction  and  leakage.  Figures  16  and  17  show  that  friction  reduced  the 
power  output  and  hence,  the  efficiency  of  the  expander  by  40%  - 50%. 
Leakage  reduced  the  efficiency  on  the  order  of  55%  - 60%.  Analytically, 
it  was  determined  that  heat  transfer  reduced  the  efficiency  approximately 
10%  - 20%. 

2)  Figures  18  through  21  show  that  the  leakage  flow  varies  only 
slightly  with  speed  but  the  ideal  flow  increases  with  speed.  Therefore, 
it  can  be  seen  from  Equation  (4.1)  that  the  reduction  in  adiabatic 
expansion  efficiency  due  to  leakage  is  decreased  at  higher  speeds. 

3)  The  frictional  power  loss  is  proportional  to  speed.  As  the 
speed  increases,  friction  tends  to  increasingly  reduce  the  efficiency. 

4)  The  frictional  power  loss  can  be  predicted  very  well,  as 
shown  in  Figure  34.  However,  the  predictions  are  dependent  on  a 
knowledge  of  the  coefficient  of  friction. 

5)  The  magnitude  of  the  leakage  flow  rates  can  be  predicted 
reasonably  well  by  the  model.  However,  prediction  of  the  correct  amount 
of  leakage  depends  on  a knowledge  of  the  leakage  paths  and  the  actual 
operating  clearances.  This  information  is  difficult  to  obtain  in 
practice. 

6)  The  accuracy  of  the  predicted  power  output  is  reduced  by 
severe  amounts  of  leakage  as  a result  of  a lack  of  information  regarding 
the  leakage  paths. 

The  purpose  of  this  study  was  not  to  develop  an  efficient  vane 
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expander,  but  to  obtain  an  understanding  of  the  phenomena  occurring 
within  the  expander.  However,  the  information  gained  provides  a basis 
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for  designing  a vane  expander  with  projected  efficiencies  50%  to  100% 
greater  than  those  measured  in  the  experimental  version. 

The  frictional  power  loss  can  be  reduced  by  decreasing  the  normal 
forces  at  the  vane-rotor  and  vane-stator  contact  points.  Supply  of  the 
steam  through  a port  in  the  end  plate  or  through  the  rotor,  at  the  vane 
base,  may  prevent  the  vane  from  being  pushed  into  the  vane  slot  by 
incoming  steam.  This  would  permit  a decrease  in  the  spring  force  at  the 
base  of  the  vanes  with  a concomitant  decrease  in  the  vane-stator  contact 
force. 

Frictional  power  losses  can  also  be  decreased  by  reducing  the 
friction  coefficient.  Injection  of  oil  around  the  periphery  of  the 
stator  would  increase  the  probability  of  obtaining  a hydrodynamic  film 
between  the  vanes  and  the  stator,  resulting  in  a decrease  in  the 
coefficient  of  friction.  An  increase  in  the  oil  flow  rate  may  result 
in  a decrease  of  the  friction  coefficient,  but  it  would  also  result  in 
more  contamination  of  the  steam.  There  may  also  be  other  material 
combinations  with  lower  coefficients  of  friction  in  the  presence  of 
steam  and  oil. 

Leakage  in  the  expander  could  be  reduced  by  matching  the 
coefficients  of  thermal  expansion  of  the  expander  components.  This 
would  result  in  lower  operating  clearances  than  were  obtained  in  the 
test  expander.  Elimination  of  the  thin  liner  which  is  susceptible  to 
distortion  could  reduce  leakage  between  the  vane  and  the  stator. 

Heat  transfer  from  the  expander  could  be  significantly  reduced 


by  elimination  of  the  pressure  transducer  and  slip  ring  cooling  adapters 
installed  in  the  test  expander.  Total  heat  losses  could  probably  be 
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APPENDIX  A 


EQUATIONS  USED  IN  THERMODYNAMIC  MODEL 


The  ordinary  differential  equations,  used  in  the  thermodynamic 
model,  for  the  pressure,  mass  and  in  the  super-heated  region,  temperature 
of  the  working  fluid  were  the  same  as  those  given  by  Wolgemuth  and  Olson 
in  Reference  1.  They  will  be  listed  here  in  dimensionless  form. 

Conservation  of  mass  for  the  control  volume  shown  in  Figure  38 
(volume  between  adjacent  vanes  plus  the  volume  of  the  pressurization 
ports  and  the  volume  of  the  trailing  vane  slot  underneath  the  vane)  shows 


M*  - m*  - m* 


(A.  1) 


where  the  subscript  (1)  Indicates  "in"  and  the  subscript  (2)  Indicates 
"out".  The  superscript  (*)  indicates  a dimensionless  variable.  The 
dimensionless  variables  are  defined  in  terms  of  the  dimensional  variables 
in  Table  3.  For  the  inlet  and  exhaust  ports 


BA* 

nlor2  “ (v/vQ) 


(A. 2) 


where  B , the  breathing  number  is  defined  as 


6°  Cp  ^ 223.8 /tT 


2TTNZ  R 


(A. 3) 
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TABLE  3 


Definitions  of  Non-Dimensional  Variables 


0 « 

D 60  c 


(Ro2Z) 


ZR  2 /v 
o o 


2TTNZR 


V ZR  2 
o 


V*(f, 


ZR  2iTNh 
o o 


P ZF 
o o 


27rRo  ZPORT 
NUMBER  OF  VANES 


The  variation  of  A with  time  is  shown  in  Figure  39  for  the  exhaust 

port,  and  for  the  inlet  port  when  the  arc  of  admission  is  greater  than 

the  arc  (angle  between  adjacent  vanes).  Figure  40  shows  the  variation 
* 

in  A for  the  inlet  port  when  the  arc  of  admission  is  less  than  arc. 
For  the  superheated  region,  the  differential  equation  for  the 


dimensionless  pressure  is 


**  •*  h **  h , ,AD.  * ,9h, 

Q + [m*  - m*  - - M -]-(_)  M (^) 


r 3h  , . n .9hx  V /. 

M 9P  t + AB%) p ’ Vo  iF  ( 


o o 

V*  ,144X14.6959, 


14.6959  h 


where 


•*  •*  . V v 

V ~ M (—*) 

o 

0.0160185  ,4.55504  , 9b' 


4.55504T  9b' 


AB  - 


4.55504  9bJ_ 

P 9T 


Based  on  the  Keenan  and  Keyes  [14]  equation  of  state  for  steam 


v - 0.0160185(4.55504  T/P  + b')  , 


the  dimensionless  temperature  in  the  superheated  region  is  computed  by 


integrating 


101 


0 XEX  XAD 

ANGULAR  POSITION  OF  LEADING  VANE 


Figure  39. 


Variation  of  Exhaust  and  Intake  Area  for  Arc  of 
Admission  and  Arc  of  Exhaust  Greater  then  Arc 
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(A.  8) 


•* 


In  the  two-phase  region,  the  P equation  is 


•* 

P 


[Q"  + m*  - m*  xr M*  ^ ~ VM*)  ] * (A*9> 


. •* 


2 h 


v,  h ' o 
fgo 


VL 

M [AC  " Vo  * 778  h 

M o 


where 


AC 


d,  , d,  h,  r-  dll,  dU  -i 

fg  l j 


10) 


l 

E1, 

»> 


Since  P and  T are  not  independent  in  the  two-phase  region,  an 
•* 

expression  for  T is  not  needed  in  that  region. 

The  differential  equation  for  the  volume  of  the  control  volume 
employed  in  Reference  1 was  modified  to  account  for  the  changing  volume 
in  the  trailing  vane  slot. 

For  the  vane  slot  and  pressurization  ports  in  Figure  41,  the 
volume  is 


V = SW  • Z • (SD  - LV  + HR)  + VH  + VOP  , (A. 11) 

bL 


where 


HR 


2 , 1/2 

-R±  + e cos(0  - ARC)  + Rq[1  - “ sin  (9  - ARC)]  .(A. 12) 

R 

o 


A 

«a 


is  the  angular  position  of  the  leading  vane,  VH  is  the  volume  of 
the  pressurization  ports  and  VOP  is  the  volume  of  the  groove  at  the 


base  of  the  vane  slot.  Since 


.» 


* 

V SL 

‘ VSL/ZRo2  • 

(A. 

then 

•* 

dV*  _ dV*  dV 

SW  dHR 

(A. 

V SL 

d0  = dV  d6 

2 d0  * 

O 

or 

•* 

V SL 

sw  r 

0 

2 

sin(0  - ARC)-  sin(0 

0 

- ARC)cos(0 

2 

- ARC) [ 1 - ~ 
R Z 

O 

sin2 (0 

-1/2-1 

- ARC)] 

(A. 

Therefore,  the  differential  equation  for  the  dimensionless  volume 
becomes 


• * 

V «*  (1  - R 


')2  cos2e  + (1  - R')  cos  e/l  - (1  - R’)2  sin2e  - (1  - R 


cos2  (8  - ARC)  - (1  - R’)  cos  (6  - ARC)/  1 - (1  - R’)2  sin2 (6 


•* 
+ V 


SL  ’ 


(a.: 


where  R'  ■ R. /R 
l o 

The  dimensionless  volume  of  the  control  volume  at  the  position 
when  9=0  (the  initial  volume)  is 
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* 

VI 


2 2 

R - R 2 R 2 

(— — ^ — — )ARC  + y-  sin (2  ARC)  + y e sin  ARC(1  - — • sin 

R 

o 


ARC)1/2  + ~y  sin'1^-  sin  ARC 

o 


+ SW(SD  - LV  + HR) 


e 


7tt 

4 


+ (OP) (OPD) 


(A. 17) 


The  differential  equation  for  the  dimensionless  power  is 


'* 

W 


* * 
P V 


(A. 18) 


In  the  superheated  region,  Equations  (A.l),  (A. 4),  (A. 8),  (A. 16) 
and  (A. 18)  where  integrated  on  the  digital  computer  using  the  IBM  360 
Scientific  Subroutine  Package  (360A-CM-03x)  version  III  subroutine  DHPCG. 
Subroutine  DHPCG  uses  the  modified  Hamming's  Predictor  Corrector  method 
[7]. 

In  the  two  phase  region.  Equations  (A.l),  (A. 9),  (A. 16),  and  (A. 18) 

were  integrated  on  the  digital  computer  using  subroutine  DHPCG. 

The  solution  of  the  equations  was  started  at  0=0  with  initial 

guesses  for  working  fluid  properties  in  the  control  volume  and  the 

* 

dimensionless  work  initialized  at  W = 0 . The  equations  were  integrated 
through  one  revolution  and  the  final  values  of  the  properties  used  as 
initial  values  for  the  next  revolution.  The  process  continued  until  the 
initial  and  final  property  values  were  the  same. 

In  the  two  phase  region,  the  steam  thermodynamic  properties  and 
their  derivatives  were  read  into  the  computer  as  tabular  data.  The 
properties  and  their  derivatives  were  obtained  from  Reference  14.  A 
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three  point  Interpolation  scheme  gave  good  values  of  the  properties  and 
their  derivatives  from  the  tabular  data. 

In  the  superheated  region,  the  steam  thermodynamic  properties  were 
obtained  from  PSU  IBM  360  computer  subroutines  based  on  References  14, 

15,  and  16.  The  partial  derivatives  of  "b"'  and  the  enthalpy  which 
appear  in  Equations  (A. 4)  through  (A. 7)  were  evaluated  from  analytic 

I 

expressions  resulting  from  differentiating  the  expressions  for  the 
enthalpy  and  "b"'  given  in  Reference  15  and  16. 
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APPENDIX  B 

DERIVATION  OF  EXPRESSIONS  FOR  TANGENTIAL  AND  NORMAL  ACCELERATIONS 

The  tangential  acceleration  of  the  vane  center  of  gravity  may  be 
written  as: 


f i 

* . * 


acg  t Ap  t + ACG/p  t + acoriolis  t 


(B.  1) 


The  magnitude  of  the  tangential  acceleration  of  a point  P(A  ) on  a 

P t 

frame  moving  with  the  rotor  at  a radius  coinciding  with  the  radius  of 
the  vane  center  of  gravity  is  in  general 


Ap  t r ROTOR  * 


(B.  2) 


However,  at  steady  state,  the  speed  of  the  rotor  is  constant  and  the 

angular  acceleration  (tt1  ) is  zero.  It  has  been  assumed  that  the 

KU1UK 

acceleration  of  the  vane  relative  to  the  moving  frame,  in  the  tangential 
direction,  ^G/p  t is  zero*  Therefore,  Equation  (B.l)  reduces  to 


Sx  t “ ACORIOLIS  t “ X VCG/P  ‘ 


(B.  3) 


Figure  3.1  indicates  that  the  distance  HR  describes  the  position  of 
the  blade  tip  relative  to  the  rotor.  From  Reference  1 


HR  = -R,  + e cos  0 + 8(1  * Ar  sin20)1/2 

1 0 R 2 

o 


(B.4) 


where  R.  is  the  radius  of  the  rotor,  R is  the  radius  of  the  stator 
i o 

and  e “ Rq  - R^  is  the  eccentricity.  Therefore, 


dHR  „ dHR  d_6 
dt  d6  dt 


(B.  5) 


27rN[-e  sin  0 - -f—  sin  0 cos  0(1  — sin20)  ] , 
o R 2 

° (B. 6) 


The  velocity  vCG/p  is  the  radial  direction.  Hence,  the  resulting 
expression  for  the  magnitude  of  the  tangential  acceleration  of  the  vane 


2 e2  e2  2 "1/2 

2(2nN)  [-e  sin  0 - — sin  0 cos  0(1 r-  sinz0)  ] 

K n 

o R 

° (B. 7) 


The  direction  of  ACG  t is  found  by  rotating  the  vector  vCG/p  90°  in 
the  sense  of  rotation  of  the  rotor. 

The  normal  acceleration  of  the  vane  is 


A__  * A + A_ _ , 

CG  n p n CG/p  n 


(B.8) 


The  magnitude  of  the  normal  acceleration  (Ap  of  a point,  on  a frame 
moving  with  the  rotor,  at  a radius  corresponding  to  the  radius  of  the 


vane  center  of  gravity  is 


(B.9) 


This  acceleration  is  always  directed  toward  the  center  of  the  rotor. 
The  magnitude  of  the  acceleration  of  the  vane  relative  to  the  moving 
frame  (\:G/p  *-s  obtained  by  differentiating  Equation  (B.6)  with 
respect  to  time.  Therefore, 


CG/p  n 


2 e2  e2  2 "1/2 

(2ttN)  i-e  cos  0 - — cos  29(1 ^ sin  0) 

o R 1 

o 


.4  2 -3/2 

j-  sin  20(1  - ~ sin  0)  ] 

R 4 R 

o o 


(B. 10) 


Vector  addition  of  the  quantities  in  Equation  (B.9)  and  (B.10)  gives  the 
normal  acceleration  of  the  vane. 
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APPENDIX  C 


DERIVATION  OF  FINITE  DIFFERENCE  EXPRESSIONS 


The  transformation  of  the  three  dimensional  Laplacian 


u2t  - 3^1  + I SI  . J_  32T  , 92T 

V 1 . 2 r 3r  + 2 ..2  + . 2 

dr  r 30  3z 


(C.l) 


into  finite  difference  form  was  accomplished  by  the  method  of  Taylor's 
Series  expansion  [11,13], 

For  interior  nodes,  as  shown  in  Figure  42,  the  following  Taylor's 
Seriea  expansions  were  written  about  the  point  i.j.k 


i-l,J,k 


i+l,J,k 


i.j+l.k 


1,3-l.k 


i.j  *k~l 


i.j  »k+l 


T _ A0  3T  . A02  a^T  _ . 

1,J,k  99  i,j,k  _5T'  302  i,j,k3!  303  i.j.k’ 


(C.2) 


Afi  9T|  A02  9^T  A0^  H T 

i.j.k  96  . . . "ZT  2 , . . 3!  -fi3  . , . * * ' 

'i.j.k  96  i.j.k  96  i,j,k 

(C.  3) 

T + Ar  9T  + Ar2  + AJL.  33T  + 

i.j.k  9r  i.j.k  9r  i.j.k 

(C.4) 

t _ Ar  H + Ar^  S^T  _ Ar^  j^T  + 

i.j.k  9r  1J>k  21  ^2  liJfk3!  9r3  1>1>k 

(C.5) 


A03  93t| 


T . . + Ar 
i.j.k  9r 


i.j.k 


x _ A. ii  + 

i.j.k  9z  i.j.k  2!  3z2  1>jk3l  ^3  i>Jtk 


(C.  6) 


- T +Az3I  + As! + As! 

i.j.k  9z  t j>k  2!  3z2  3J  9j,3  ± ^ ^ >R 


+ 


(C.7) 
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It  is  necessary  to  add  Equations  (C.2)  through  (C.7)  in  a way  that 
will  permit  the  attainment  of  a finite  difference  approximation  to 
Equation  (C.l).  Therefore,  Equations  (C.2)  through  (C.7)  were  multiplied 
by  the  coefficients  A,B,C,D,E,  and  F , respectively.  The  following 
constraints  were  applied  to  determine  the  coefficients: 


Z coefficients  of  — = 0 


Z coefficients  of 


Z coefficients  of 


Z coefficients  of  -P-  = — 
9r  r 


Z coefficients  of 


Z coefficients  of 


(C.8) 


The  coefficients  determined  by  application  of  those  constraints  are: 


l/(rA9)‘ 


2r  + Ar 

2rAr2 

2r  - Ar 
2 


- F - 1/Az' 


(C.9) 


I 


Multiplication  of  Equations  (C.2)  through  (C.7)  by  the  constants  A 
through  F , respectively,  permits  the  attainment  of  the  following 
approximation  for  Equation  (C.l): 


y^T  _ i n + ^r—  — f 1 X 

(rA0) 2 i_1,j,k  2rAr2  1,:}’k  (rA0)2  1+1,j,k 


+ .^r  _ Ar  j + JL  f + _1_  t 

a 2 i.j.k+1  A 2 Ai,J,k-l 


2 2r  + Ar  2r  - Ar  2 
(rA0)2  2rAr2  2rAr2  Az2  i,j,k 


For  grid  points  on  a boundary  as  shown  in  Figure  43,  where  heat 

generation  and  convection  are  present,  Equation  (C.l)  was  approximated 

by  using  the  fictitious  node  T 1 , . The  boundary  condition  on  the 

— 

surface  was r — = h(T  - T ) + q"  . (C.ll) 

3r  °°  w Hgen 

The  derivative  on  the  left  hand  side  of  Equation  (C.ll)  was 
approximated  using  a central  difference  scheme.  Therefore, 


Ti,j-l,k  _ Ti,_1+l,k  _ h _ ^gen 

2Ar  k k °°  i,j  ,k'  k 


_ 2Arh,  T n , I..  2Ar 

i.J-l.k  k u°°-ii,j,k;  qgen  k i,j+l,k 


substitution  of  Equation  (C.12)  into  Equation  (C.10)  gives 
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2 1 
V x = ±—7T  T 


+ ^4lT. 


T. 


(rA6)2  1 1*j>k  2rAr2  i>J+1>k  (rA6)2  1+1*j’k 


Az 


2 "i.j.k-l  Az2  i,j,k+l 


+ 2r  - Ar  T 


2rAr 


2 i,j+l,k 


+ (2r-~  Ar)  JL  t + 2r  - -Ar  Vn 

k 2rAr  ' k « 2rAr  k 


2 2r  + Ar  2r  - Ar  2 2r  - Ar  h 1 
, 2 + o a . 2 o,2  + ,„2  + 1 rAr  J k 


[ (rA0)  2rAr  2rAr  Az 


i.j.k 
(C. 13) 


i' 

f 

; 


n 


r 


The  fictitious  node  technique  as  described  above  was  applied  to  all  the 
boundaries. 

In  the  steady  state  analysis,  Equation  (C.l)  is  equal  to  zero. 
Therefore,  the  temperature  at  each  node  can  be  determined  by  using 
Equations  like  (C.10)  and  (C.13).  For  example,  using  Equation  (C.10) 
the  expression  for  T . is 


i.J.k 


+ 2r  T 


<rA6) 


2 *i-l,j,k  2rAr2  i,j+l,k  (rAQ)2  i+l,j,k 


+ ^#T 


2rAr2  Az2  Az2 


2 2r  + Ar  + 2r  - Ar  + _2_ 
[ (rA9)2  2rAr2  2rAr2  Az2 


(C. 14) 


Solving  for  T at  each  node,  results  in  a set  of  algebraic 

1 » J » ^ 

expressions  for  T which  must  be  solved  simultaneously.  This  was 

i « J • k 

accomplished  using  the  Gauss-Seidel  iterative  method  with  the  incorpora- 
tion of  Young's  over  relaxation  factor  [7]. 


J 
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For  the  transient  analysis,  the  equation 


V2T  = I 31 

a 3t 


(C. 15) 


was  approximated  by  finite  differences.  For  the  explicit  solution,  a 

bt 

finite  difference  expression  for  was  obtained  using  a forward 

difference.  Therefore,  the  approximation  was 


tN+1  - tN 
3T  _ i,1,k  ~ i,i,k 

3t  At 


(C . 16) 


h i 


The  superscripts  (N+l)  and  (N)  refer  to  time  (t  + At)  and  time  (t), 
respectively. 

The  finite  difference  approximation  of  Equation  (C.15),  for  an 
interior  node  is: 


TN+1  n 
thlA  ~ i|J,k 

aAt 


+ 2r  + Ar  TN 


(rA0) 


2 i-1, j ,k  2rAr2  i.j »k+l 


* I 


+ -4tn 


(rA9) 


2 i+l,j,k 


+ 2r_^r  ^ 


+ -4  tn 


+ -4tn 


2rAr 


2 *i,j-l,k  Az2  *i,j,k+l  a,2  ‘i.j.k-l 


Az 


2 2r  + Ar  + 2r  - Ar  _2_ 
{ (rA9)2  2rAr2  2rAr2  Az2 


rN 

l,j,k 


(C.17) 


The  approximations  to  Equation  (C.15)  for  each  node  can  be  solved  for 

are  all  known,  iteration  is  not 
required  to  solve  the  set  of  algebraic  equations. 


N+l  N 

T.  , . Since  the  values  of  T,  . . 
i,j,k  1 * j >k 
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APPENDIX  D 
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DERIVATION  OF  FRICTIONAL  HEAT  FLUX  EQUATIONS 

The  total  rate  of  frictional  heat  generated  at  the  interface  of 
two  semi-infinite  solids  in  sliding  contact  is 


q 


II 

TOTAL 


qV  + q' 


(D.l) 


where  q'^  and  are  the  heat  transfered  into  body  I and  II,  respectively. 

The  heat  conduction  equation  for  each  body  is 


32T  1 3T 

9z2  = a 3t 


(D.  2) 


The  boundary  conditions  for  each  body  are: 


BODY  I 


z = o 


3T, 


1 9z 


q 


n 

l 


where 
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BODY  II 


z = 0 


, 9T  _ ■„ 

k2  3z  q2 


z = -oo 


t2  i 


(D.4) 


The  initial  condition  for  each  body  is  at 


t = 0 


T = 0 . 


(D.5) 


Taking  the  Laplace  transform  of  Equation  (D.2)  gives 


rl-S  t<x.«)  - o • 


(D.b) 


The  boundary  condition  at  z = 0 for  solid  I becomes 


Oi  / \ * 

■^(°’s)  = k^ 


(D.  7) 


when  is  assumed  to  be  constant.  The  boundary  condition  at  z = 00 

becomes 


T(°°,s)  + » . 


(D.8) 


The  general  solution  to  Equation  (D.6)  is 


s/a  z 


T(x, s)  » C^e  + C^e 


(D.9) 
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Applying  boundary  conditions  (D.7)  and  (D.8)  gives 


T(z,s) 


/ A 11  ^ ^ / Ct  c 

,ql,2rt  -3/2  /Ctl 

(k[)  ai  s 


From  Reference  17  the  inverse  transform  of  Equation  (D.10)  is 


- (z  /a  ) /At  /T 


/M,  7 / A. 

<ir>  \ Vi  e - vV  >erfc 


Similarly  for  body  II 


T(z,0)  = 


2 

(,2)  a2 

k 


o t 

Vi  e 


■U'/apMt  /% 

- ( ) erfc  


Assuming  no  contact  resistance 


7^(0, t)  = T2(0,t)  . 


Therefore,  by  equating  Equations  (D.ll)  and  (D.12),  it  is  found  that 


h a_i 

ql  q2  k j/  ' 


Combining  Equations  (D.l)  and  (D.14)  gives 


TOTAL 

Pp 

V “i 


mtotal 
k m - 
1 + rV  — 

ky  a. 
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